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SUMMARY 


An  analytical  method  for  predicting  the  stability  characteristics  of  tilt¬ 
wing  VTOL  aircraft  in  the  transition  speed  range  is  presented.  Sample 
calculations  based  on  an  assumed  tilt-wing  VTOL  transport  configuration  of 
the  XC-142A  class  with  double  slotted  flaps  are  given. 

Particular  emphasis  is  placed  on  the  sensitivity  of  the  results  to  various 
assumptions  made  in  the  analysis.  The  contributions  of  the  various  air¬ 
craft  components  and  the  aerodynamic  interactions  of  the  components  to  the 
stability  derivatives  are  discussed,  as  well  as  the  changes  in  the  charac¬ 
teristic  modes  of  motion  of  the  vehicle  that  result  from  variations  in  the 
stability  derivatives. 

The  trim  conditions  of  the  vehicle  are  shown  to  be  quite  sensitive  to  the 
prediction  of  the  flap  characteristics. 

A  limited  comparison  of  the  calculated  results  with  experimental  data  ob¬ 
tained  from  a  dynamic  model  of  the  XC-142A,  which  is  somewhat  dissimilar 
from  the  assumed  configuration,  is  presented,  lhis  comparison  indicates 
that  the  trends  of  the  stability  derivatives  are  correctly  predicted.  Hie 
agreement  between  theory  and  experiment  is  good  in  hovering;  however,  as 
the  wing  incidence  is  reduced,  the  difference  between  theory  and  experi¬ 
ment  becomes  quite  large. 


FOREWORD 


This  research  was  performed  hj  the  Department  of  Aerospace  and  Mechanical 
Sciences,  Princeton  University,  under  the  sponsorship  of  the  United  States 
Army  Aviation  Materiel  Laboratories  Contract  DA  44-177 -AMC-8(T),  with 
financial  support  from  the  United  States  Navy,  Bureau  of  Weapons,  and  thv. 
United  States  Air  Force  Flight  Dynamics  Laboratory  of  the  Research  and 
Technology  Division.  The  research  was  monitored  by  Mr.  William  E.  Sickles 
of  the  United  States  Army  Aviation  Materiel  Laboratories. 

Hie  research  was  conducted  by  G.  Beppu  and  Assistant  Professor  H.  C. 
Curtiss,  Jr.,  of  Princeton  University. 


v 


TABLE  OF  CONTENTS 


SUMMARY 
FOREWORD 
LIST  OF  FIGURES 
LIST  OF  SYMBOLS 
INTRODUCTION 

METHOD  OF  ANALYSIS  AND  ASSUMPTIONS 

DISCUSSION  OF  RESULTS 

CONCLUSIONS 

RECOMMENDATIONS 

REFERENCES 

DISTRIBUTION 

APPENDIXES 

A.  Physical  Parameters  of  Assumed  Aircraft  and 
Description  of  Apparatus  and  Experiments 

B.  Analytical  Expressions  for  Stability  Derivatives 


vii 


s 


Page 

iii 

v 

viii 

xi 

1 

2 

10 

32 

34 

66 

68 


69 

76 


LIST  OF  FIGURES 


1  Trim  Velocity  vs.  Wing  Incidence  Angle  for  Various 
Values  of  Flap  Effectiveness 

2  Propeller  Hi  rust  and  Wing  Lift  vs.  Wing  Incidence 

3  Propeller  Blade  Angle  and  Nondimensional  Induced 
Velocity  vs.  Wing  Incidence 

4  Slipstream  Velroity  vs.  Wing  Incidence  for  Various 
Values  of  Flr.p  Effectiveness 

5  Effective  Wing  Angle  of  Attack  vs.  Wing  Incidence 

6  Pitching  Moment  Required  TO  Trim  vs.  Wing  Incidence 

7  Induced  Velocity  Derivatives  vs.  Wing  Incidence 

8  Velocities  at  Propeller  Plane  vs.  Wing  Incidence 
as  a  Function  of  Trim  Velocity 

9  Dependence  of  Various  Slipstream  Parameters  on  Wing 
Incidence/Trim  Velocity  Relationship 

10  Rate  of  Change  of  Slipstream  Velocity  With  Horizontal 
Velocity  vs .  Wing  Incidence 

11  Rate  of  Change  of  Effective  Wing  Angle  of  Attack  With 
Horizontal  Velocity  vs.  Wing  Incidence 

12  Rate  of  Change  of  Thrust  With  Horizontal  Velocity  vs. 
Wing  Incidence 

13  Rate  of  Change  of  Wing  Lift  With  Horizontal  Velocity 
vs.  Wing  Incidence 

14  Stability  Derivatives  -  Rate  of  Change  of  Horizontal 
Force  With  Horizontal  Velocity  vs.  Wing  Incidence 


Page 

35 

36 


37 


38 

39 

40 

41 

42 


*6 

44 


15  Stability  Derivatives  -  Rate  of  Change  of  Vertical 

Force  With  Horizontal  Velocity  vs.  Wing  Incidence  45 

16  Rate  of  Change  of  Slipstream  Velocity  With  Vertical 

Velocity  vs  •  Wing  Incidence  46 


vlii 


V. 


* 

~  V 

f 

"\ 

’  *  '■  r * 

f 

Figure 

Page 

17 

Rate  of  Change  of  Effective  Wing  Angle  of  Attack 

With  Vertical  Velocity  vs.  Wing  Incidence 

46 

18 

Rate  of  Change  of  Propeller  Thrust  With  Vertical 
Velocity  vs.  Wing  Incidence 

19 

Rate  of  Change  of  Wing  Lift  With  Vertical  Velocity 
vs.  Wing  Incidence 

47 

20 

Stability  Derivatives  -  Rate  of  Change  of  Horizontal 
Force  With  Vertical  Velocity  vs.  Wing  Incidence 

48 

21 

Stability  Derivatives  -  Rate  of  Change  of  Vertical 
Force  With  Vertical  Velocity  vs.  Wing  Incidence 

49 

22 

Stability  Derivatives  -  Rate  of  Change  of  Pitching 
Moment  With  Horizontal  Velocity  vs.  Wing  Incidence 
(No  Horizontal  Tail) 

50 

23 

Horizontal  Tall  Contributions  to  Pitching  Moment 

Change  With  Horizontal  Velocity 

51 

24a 

Stability  Derivatives  -  Rate  of  Change  of  Pitching 
Moment  With  Vertical  Velocity  vs.  Wing  Incidence 

52 

24b 

Wing  Contributions  to  Rate  of  Change  of  Pitching 

Moment  With  Vertical  Velocity  vs.  Wing  Incidence 

53 

25 

Rate  of  Change  of  Pitching  Moment  With  Pitching 
Velocity  vs  •  Wing  Incidence 

54 

26 

Locus  of  Roots  as  a  Function  of  Wing  Incidence 

55 

27 

Locus  of  Loots  as  a  Function  of  Wing  Incidence  - 
Sensitivity  to  Xy 

56 

28 

Locus  of  Roots  as  a  Function  of  Wing  Incidence  - 
Sensitivity  to  Zy 

57 

29 

Locus  of  Roots  as  a  Function  of  Wing  Incidence  - 
Sensitivity  to  Mg 

58 

30 

Locus  of  Roots  as  a  Function  of  Wing  Incidence  - 
Sensitivity  to  My 

59 

31 

Locus  of  Roots  as  a  Function  of  Wing  Incidence  - 
Sensitivity  to  My 

60 

i 


Figure 

Page 

32 

Comparison  of  Theory  and  Experiment 
vs.  Wing  Incidence 

-  Trim  Velocity 

6l 

33 

Comparison  of  Theory  and  Experiment 
Incidence 

-  Xu  vs.  Wing 

34 

Comparison  of  Theory  and  Experiment 
Incidence 

-  Zy  vs.  Wing 

62 

35 

Comparison  of  Theory  and  Experiment 
Incidence 

-  Mjj  vs.  Wing 

36 

Comparison  of  Theory  and  Experiment 
Incidence 

-  M»  vs.  Wing 

63 

37 

Assumed  Flap  Characteristics 

38 

Horizontal  Force  vs.  Angle  of  Attack  for  Various 

Wing  Incidences 

64 

39 

Pitching  Moment  vs.  Angle  of  Attack 
Wing  Incidences 

for  Various 

65 

A-l 

Axis  System  and  Notation 

73 

A-2 

Flap  Program  vs.  Wing  Incidence 

74 

A-3 

Aircraft  Geometric  Characte_  i sties 

A-4 

General  Arrangement  1/10  Scale  XC-142  Model 

75 

X 


,44  I 


■ft  "  5M 


A 


7 


t *  i 


'•’HWl.'lWl* 


Tty*******'  •*« 1  «aWW  WWPWft 


LIST  OF  SYMBOIfi 


Aq. 

AR 

a 

a.c. 

b 

c 

c.g. 

% 

CD,Sc 

Cto>S 


"TjS 


Propeller  disc  area  (189  square  feet) 

Tall  rotor  disc  area  (50.2  square  feet) 

Wing  aspect  ratio  (8.53) 

Propeller  blade  lift  curve  slope  (5*73  per  radian) 

Aerodynamic  center 
Number  of  blades  (4) 

Wing  span  (67.5  feet) 

Blade  chord  (feet) 

Center  of  gravity 

Drag  coefficient  at  zero  lift 

Drag  coefficient  at  zero  lift,  based  on  the  slipstream  velocity 


Wing  lift  curve  slope  in  slipstream,  based  on  slipstream  ve¬ 
locity  and  effective  angle  of  attack  (per  radian) 

Wing  lift  curve  slope  in  a  uniform  free  stream  (per  radian) 
Lift  curve  slope  of  horizontal  tail  (per  radian) 


Lift  coefficient  of  horizontal  tail 


Thrust  coefficient, 


p  tt  R2  (OR)2 


T/A 


Thrust  coefficient,  C_  e  -  - 

T,S  q  +  t/A 

Fuselage  drag  parallel  to  free  stream  (pounds) 

Wing  drag,  parallel  to  slipstream  velocity  (pounds) 

xi 


f 

H 

I 


K 


K, 


“s 

Lp 


Drag  of  horizontal  tail  (pounds) 

Equivalent  flat  plate  area  of  fuselage  (5.0  square  feet) 

Inplane  force  of  the  propeller  (pounds) 

Moment  of  inertia  of  the  airplane  about  Y-axis  (100,000 
slug-feet  squared) 

Horizontal  tail  incidence  angle,  measured  from  fuselage 
reference  line  (radians) 

Wing  incidence,  measured  from  fuselage  reference  line  (radians 
or  degrees) 

Correction  factor  for  wing  lift  curve  slope  in  slipstream  (0.5) 
Fuselage  moment  factor  (0.83) 

Wing  lift,  perpendicular  to  slipstream  velocity,  V  (pounds) 

R 

Horizontal  tail  lift  (pounds) 


'-'H’^'HT’ 

TR'*0P* 

lP’l0U,lU 

M 

Mf 

M!D 

mfd 


Various  geometric  distances  from  center  of  gravity 
(Figure  A-lb)  (feet) 


Pitching  moment,  nose  up  positive  (foot-pounds) 

Fuselage  pitching  moment  (foot-pounds) 

Dovnwash  reduction  factor,  nondimensional  (Equation  16) 

Propeller  pitching  moment  due  to  inflow  distribution  arising 
from  a  pitching  velocity  (foot-pounds) 

Propeller  pitching  moment  due  to  induced  velocity  distribution 
across  propeller  disc  (foot-pounds) 

Stability  derivatives 


=  1 
I  fcu 


Mu  -  -  -  , 


«w 


=  1  &m 

I  dw 


”9 


_  1  fcM 

’  i  ae 


m 


Mass  of  airplane  (slugs) 


xii 


'»7 

/ 


r. 


N 

q 

r 

R 

*T 

S 

S 


u 

V 

VF 

v 

V„ 


W 


f 

*u*r 


Number  of  propellers  (4) 

Free  stream  dynamic  pressure  (pounds  per  square  foot) 

Radius  of  blade  element  from  axis  of  rotation  (feet) 

Propeller  radius  (7*75  feet) 

Tail  rotor  radius  (4  feet) 

Wing  area  (53^  square  feet) 

Horizontal  tail  area  (140  square  feet) 

Propeller  thrust  (pounds) 

Tail  rotor  thrust  (pounds) 

Horizontal  velocity  perturbation  (feet  per  second) 

Flight  velocity  (feet  per  second) 

Slipstream  velocity  (feet  per  second) 

Propeller  induced  velocity  (feet  per  second) 

Average  value  of  propeller  Induced  velocity  (feet  per  second) 

Slope  of  propeller  induced  velocity  with  nondimens ional  radius 
(Equation  3)  (feet  per  second) 

Tail  rotor  induced  velocity  (feet  per  second) 

Gross  weight  of  aircraft  (37/350  pounds) 

Vertical  velocity,  downward  motion  of  the  aircraft  is  positive 
(feet,  per  second) 

Longitudinal  force,  positive  forward  (pounds) 

Longitudinal  force  due  to  fuselage,  positive  forward  (pounds) 

Stability  derivatives 


i  ax 

m  du 


i  ax 

m  dv 


xiii 


\ 


> 


i'%* 


V' 


v*. 


i 


i*^rf  > 


z  ,z 

u'  w 


cr 

or- 


0 

0. 


0.75R 


O-'ORrp 


Vertical  force,  positive  downward  (pounds) 
Stability  derivatives 


z  =  , 

U  D  3u 


Z  .1*5 


w  m  9w 

Angli  of  attack  of  the  fuselage  (radians  or  degrees) 

Angle  of  attack  for  zero  lift  of  the  wing  (radians  or  degrees) 

Effective  wing  angle  of  attack  in  slipstream  (Equation  11) 

( radians ) 

Flap  deflection  angle  (radians  or  degrees) 

Blade  section  profile  drag  coefficient  of  the  propeller  (0.016) 
Blade  section  profile  drag  coefficient  of  the  tail  rotor  (0.016) 

Down wash  angle  (radians) 

Propeller  downwash  angle  at  tail,  ep=(iw  +  a-ae)  (radians) 

Downwash  angle  at  the  tail  rotor  due  to  horizontal  tail  lift 
(radians) 

Aircraft  pitch  angle,  positive  nose  up  (radians) 

Propeller  blade  pitch  angle  at  0.75R  (radians) 

Tail  rotor  blade  pitch  angle  at  0.75R  (radians) 


Nondimens ional  induced  velocity  of  the  propeller  X^  =  ^ 
(always  positive) 


Nondimens  ional  induced  velocity  of  the  tail  rotor  X^  = 

(always  positive)  T  Vin/T 

Nondimensional  velocity  parallel  to  the  propeller  plane 


V>  = 


V  sin(iv  +  or) 

"niT 


or,  using  Equation  11,  p,  = 


VR  sin  °e 
OR 


Nondimensional  velocity  parallel  to  the  plane  of  the  tail  rotor 
V  cos  or 

(ORL 


*r  = 


xiv 


*** v«  \*ni  ,( 


V 


V 


T 


a 


% 


Y 


(  ) 

(  ) 

O 

(  ) 


HT 

T 


o 


Nondimenslonal  velocity  perpendicular  to  the  propeller  plane 


v  = 


V  COs(iw  +  cy) 

nR 


Nondimenslonal  velocity  perpendicular  to  the  plane  of  the  tail 

V  sin  0/ 
rotor  v  =  — — - — 

T  (ta)m 


Solidity  of  main  propeller  o  =  —  ,  (0.153) 

ttR 

Solidity  of  the  tail  rotor  (0.124) 


Air  density  (slugs  per  cubic  foot) 

Rotational  velocity  of  the  propeller  (129  radians  per  second) 
Rotational  velocity  of  the  tail  rotor  (363.4  radians  per  second) 


=  V(X1  +  v) 


2  2 

+  ii  nondimenslonal  velocity 


Propeller  blade  azimuth  angle  (measured  from 

Tail  efficiency  factor  T  —  ^ 

9 


at  propeller  plane 
downwind  position) 


ADDITIONAL  NOTATION 


Horizontal  tail 
Tail  rotor 

Derivative  with  respect  to  time 
Initial  value 


xv 


BLANK  PAGE 


INTRODUCTION 


Tli5  analytical  prediction  of  the  stability  and  control  characteristics  of 
tilt-wing  VTOL  aircraft  at  high  wing  incidences  with  correspondingly  low 
forward  speeds  presents  a  number  of  difficulties.  The  first  is  the  alge¬ 
braic  complexity  of  the  expressions  obtained  for  the  stability  derivatives, 
arising  from  the  numerous  components  and  aerodynamic  interactions  that  con¬ 
tribute  to  the  aerodynamic  forces  acting  on  the  vehicle.  A  second  and  more 
important  difficulty  is  associated  with  the  lack  of  suitable  aerodynamic 
data  available  on  the  various  components  for  stability  and  control  analy¬ 
sis.  In  addition,  the  nature  of  the  flow  fields  involved  makes  theoretical 
results  difficult  to  apply  since  simple  closed-form  expressions  are  gener¬ 
ally  not  available. 

The  algebraic  complexity  may  be  circumvented  by  the  use  of  high-speed 
digital  computing  equipment.  However,  the  importance  of  the  sensitivity  of 
the  results  of  an  analysis  to  various  assumptions  made  during  the  course  of 
the  analysis  becomes  clouded  by  the  length  of  the  analytical  expressions 
used. 


The  objective  of  this  study  is  to  clarify  the  aerodynamic  interactions 
entering  into  the  stability  calculations  and  to  determine  their  signifi¬ 
cance  to  the  prediction  of  the  stability  characteristics  of  the  aircraft. 


The  approach  presented  here  is  analytical.  A  simplified  mathematical  model 
of  the  propeller  and  the  wing-slipstream  interaction  is  used  to  predict  the 
trim  conditions  and  the  stability  derivatives  of  an  assumed  tilt-wing  VTOL. 
The  results  are  compared  where  possible  with  experimental  data  taken  on  a 
0.10  scale  dynamic  model  of  the  XC-142A.  The  experiments  were  conducted 
on  the  Princeton  Dynamic  Model  Track.  The  test  set-up  and  the  model  are 
described  in  Appendix  A. 

The  major  contributions  to  each  stability  derivative  are  analyzed,  and  the 
significant  effects  are  discussed  for  level  flight  trim  conditions. 

The  report  consists  of  the  following  parts.  In  the  section  entitled  Method 
of  Analysis  and  Assumptions,  the  mathematical  models  used  for  the  various 
components  of  the  vehicle  and  the  assumptions  involved  are  described.  In 
the  next  section,  the  trim  calculations  are  discussed,  and  in  particular 
the  importance  of  accurately  predicting  the  aerodynamic  characteristics  of 
flaps  is  discussed.  The  next  section  discusses  the  stability  derivatives 
themselves  and  the  important  contributing  factors  to  the  stability  deriva¬ 
tives.  Then,  the  sensitivity  of  the  modes  of  motion  of  the  vehicle  to 
variations  in  the  stability  derivatives  is  discussed.  Finally,  a  com¬ 
parison  of  the  analytical  results  with  experimental  data  is  discussed. 


1 


METHOD  OF  ANALYSIS  AND  ASSUMPTIONS 


Other  methods  or  approach  to  tilt-wing  stability  and  control  are  presented 
in  References  1,  2  and  3-  References  1  and  2  are  based  primarily  on  ex¬ 
perimental  data,  and  only  certain  aspects  of  the  problem  are  treated  ana¬ 
lytically.  The  analysis  presented  in  Reference  3  is  quite  similar  to  that 
presented  in  this  report  but  with  two  exceptions.  The  airplane  configu¬ 
ration  investigated  in  Reference  3  had  flapping  propellers,  while  in  this 
investigation  a  conventional  (rigid)  propeller  is  considered.  In  addition, 
in  Reference  3,  a  rather  extreme  analytical  model  was  taken  for  the  wing 
characteristics  to  include  the  effects  of  wing  stall.  In  the  following, 
it  is  assumed  that  the  use  of  flaps  prevents  wing  stall  and  no  variation 
in  wing  lift  curve  slope  with  effective  wing  angle  of  attack  is  included. 
This  should  not  be  an  important  limitation  on  the  analysis  since  large 
areas  of  separated  flow  on  the  wing  must  be  avoided  in  this  type  of  air¬ 
craft  in  order  to  obtain  satisfactory  flying  qualities  (Reference  9). 

The  sample  calculations  in  this  report  are  based  on  a  typical  four-pro¬ 
peller  tilt-wing  VTOL  configuration.  The  physical  parameters  of  the  ve¬ 
hicle  are  given  in  Appendix  A.  Simplifying  assumptions  are  used  wherever 
possible,  and  are  discussed  in  detail  in  the  following  sections. 

PROPELLER  FORCES  AND  MOMENTS 


The  propeller  forces  and  moments  are  treated  by  methods  that  are  similar 
to  those  used  for  the  helicopter  rotor  in  stability  and  control  investi¬ 
gations  (Reference  11,  chapter  20,  and  Reference  13).  The  primary 
assumptions  are: 

1.  The  average  value  of  the  induced  velocity  at  the  propeller 
plane  is  given  by  momentum  theory  (Reference  10,  page  185) 
and  is  expressed  as  follows: 


o 

QR 


2  +  v)  +  p," 


2? 


(1) 


The  mass  flow  influenced  by  the  propeller  is  based  on  the 
resultant  velocity  at  the  rotcr  plane  (Reference  10).  The 
inplane  component  of  the  free -stream  velocity,  y.,  is  re¬ 
tained  in  this  analysis  rather  than  neglected  as  in  Refer¬ 
ences  2  and  '[ .  The  validity  of  neglecting  this  component  of 
velocity  at  the  rotor  plane  depends  upon  the  "wing  incidence/ 
forward  speed"  relationship  (References  1  and  2).  Neglect 
of  this  component  will  significantly  influence  the  propeller 
derivatives  at  high  wing  incidences,  if  the  relative  forward 
speed  is  high. 


2 


2.  The  propeller  induced  velocity  is  assumed  to  be  uni  Torn,  in  a 
plane  perpendicular  to  the  direction  of  motion  and  to  in¬ 
crease  linearly  1'rom  the  leaning  edge  ol“  the  propeller  disc 
to  the  trailing  edge  in  the  rash  ion  given  ly  Reference  8, 
that  is , 

v  =  vu  +  V1  (jp  cos  Y  (r-) 

where 

V  s  i  n("  i  +  a) 

v  =  v  — — — — — — —  ^j) 

1  0  V  cos ( i  +  o)  +  v 

W  o 

and  v  is  obtained  from  Equation  (l). 

Wliilo  this  is  a  rough  approximation  to  the  actual  induced 
velocity  variation  across  the  propeller  disc,  this  variation 
is  the  source  of  the  propeller  pitching  moments.  Unil'om 
induced  velocity  will  result  in  zero  propeller pitching 
moments  in  conflict  with  experimental  data.  See,  for  ex¬ 
ample,  Reference  7- 

3.  The  propeller  blades  are  assumed  to  be  infinitely  rigid  in 
bending  and  torsion. 

4.  It  is  assumed  that  the  effects  of  twist  and  taper  of  the 
blades  may  be  taken  into  account  by  use  of  an  equivalent 
pitch  angle  and  weighted  solidity  (Reference  10,  page  86). 

The  blade  pitch  angle  used  in  the  following  formulae  is 
taken  to  be  that  at  75-percent  radius. 

5.  The  influence  of  the  presence  of  the  wing  on  the  propellers 
and  mutual  interference  between  propellers  is  neglected. 
Limited  experimental  data  indicate  that  the  presence  of  the 
wing  may  have  a  strong  influence  on  the  propeller  (Reference 
7);  however,  there  is  not  sufficient  information  to  estimate 
the  importance  of  this  effect.  It  would  be  expected  that 
interference  between  propellers  may  be  important  in  computing 
lateral  stability  derivatives,  but  not  in  longitudinal  calcu¬ 
lations  . 

While  the  above  assumptions  yield  a  highly  simplified  model  of  the  actual 
propeller  and  may  not  be  satisfactory  for  performance  estimation  or  compu¬ 
tation  of  propeller  blade  stresses,  unpublished  calculations  by  one  of  the 
authors  indicate  that  the  basic  trends  of  propeller  aerodynamics  are  ade¬ 
quately  represented  for  stability  and  control  investigations. 

On  the  basis  of  these  assumptions,  blade  element  equations  may  be  developed 
to  evaluate  propeller  thrust,  inplane  force,  and  pitching  moment.  The 
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results  of  Reference  10,  pages  190  and  1)8,  may  be  used  lor  calculating 
thrust  and  inplane  force.  The  flapping  terms  in  the  inplane  force  ex¬ 
pression  are  set  equal  to  zero. 


?  .  ^ 2  aa  reO-75R  .  3  Jx  Xt  +  v1 
T  =  p  n  R  (QR)  —  I  ^  (1  +  -  p,  )  -  2 


o  o  rt  (X ,  +  v)  “1 

H  =  p  n  R  (QR)  j-  +  :  eo.75Rj  ^ 


(b) 

(5) 


The  expression  for  propeller  pitching  moment  due  to  the  variation  in 
Lniuced  velocity  across  the  disc  may  be  obtained  from  Reference  If  as 


(M>ro  =  m  11  n  ° a  r3  vi  (6) 

To  compute  the  propeller  pitching  moment  due  to  a  pitching  velocity  about 
the  propeller  hub,  the  change  in  inflow  distribution  arising  from  this  ve¬ 
locity  is  taken  into  account.  The  induced  velocity  is  assumed  to  be  inde¬ 
pendent  of  pitching  velocity.  Physically,  it  is  clear  that  there  would  1  e 
a  redistribution  of  induced  velocity  due  to  the  change  in  load  distribution 
as  a  result  of  angular  motion.  This  induced  velocity  variation  would  act 
to  re  iu  e  the  rotor  damping  as  compared  to  the  value  bared  on  constant 
induced  velocity.  Therefore,  a  somewhat  larger  value  of  damping  will  be 
obtained  than  would  be  present.  However,  as  seen  in  the  following,  the 
propeller  angular  damping  contribution  is  small  and  so  the  significance  of 
-my  errors  incurred  would  be  small.  Then  the  theoretical  expression  for 
the  contribution  of  the  propeller  to  the  pitch  damping  becomes 


r  =  _  £_S2_U  QR5  (7) 

v  de  'ID  16 

For  pitching  moment  trim  at  low  speeds,  the  examnle  airplane  has  a  tail 
rotor.  The  contributions  to  the  stability  derivatives  Iran  the  tail  ro  or 
are  computed  usirg  the  same  assumptions  as  for  the  propellers.  The  down- 
\.u.d.  from  the  wing  and  horizontal  tail  is  included  in  computation  oi  the 
inflo  velocity,  XT>  of  the  tail  rotor. 

XT  =  (7^  [-  V  sin  <y  +  eV  +  ^  V  +  (8) 

where 

2Cl 

*T  =  TTa^  (or  +  %  ■  ^ 

The  inplane  force  and  pitching  moment  of  the  tail  rotor  are  neglected. 
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WING  LIFT  AND  DRAG 


The  wing  lift  and  drag  are  computed  taking  into  account  what  are  considered 
to  be  the  first-order  effects  of  the  presence  of  the  propeller  slipstream. 
On  the  example  airplane,  the  wing  is  entirely  immersed  in  the  slipstream. 
The  following  assumptions  regarding  the  nature  and  influence  of  the  slip¬ 
stream  are  made: 

1.  The  slipstream  velocity  is  computed  on  the  basis  of  a  uniform 
propeller  induced  velocity  far  downstream  2vQ.  This  is 

considered  reasonable  since  the  contraction  of  the  lipstream 
takes  place  a  short  distance  downstream  of  the  disc  (Refer¬ 
ence  19) •  However,  certain  of  the  derivatives  computed  in 
the  following  are  sensitive  to  this  assumption,  and  it  may  be 
desirable  to  use  a  value  somewhat  less  than  2  to  account  for 
the  contraction  and  downstream  dissipation  of  the  slipstream. 

2.  The  effects  of  slipstream  rotation  are  neglected.  While 
slipstream  rotation  may  have  a  significant  effect  on  spanwise 
load  distribution  on  the  wing  (Reference  1),  it  is  assumed 
that  its  influence  on  the  longitudinal  stability  character¬ 
istics  will  not  be  important  unless  large  areas  of  separated 
flow  result. 

This  assumption  is  implicit  in  the  usual  method  of  data 
presentation  where  it  is  assumed  that  the  aerodynamic  forces 
on  a  wing-propeller  combination  are  only  a  function  of  a 
and  C,p  s  (References  7  and  18  for  example),  Hiat  is,  the 

slipstream  rotation  depends  upon  the  torque  of  the  propeller 
and  thus  both  on  advance  ratio  and  blade  angle.  Hius,  there 
is  li*-+le  or  no  aerodynamic  data  to  verify  the  validity  of 
this  asjumption.  Experiments  should  be  conducted  at  differ¬ 
ent  combinations  of  blade  angle  and  advance  ratio  that 
produce  the  same  propeller  thrust  at  different  levels  of 
torque  to  determine  the  importance  of  slipstream  rotation. 

3.  The  effects  of  the  finite  extent  of  the  slipstream  flow  are 
taken  into  account  by  a  reduction  in  the  wing  lift  curve 
slope.  It  is  shown  in  Reference  15  that  the  finite  extent  of 
the  slipstream  will  reduce  the  lift  curve  slope  of  the  wing 
based  on  slipstream  velocity,  compared  to  its  value  in  a 
stream  of  infinite  extent  of  the  same  velocity. 

A  simplified  approach  such  as  presented,  for  example,  in 
Reference  5j  is  considered  satisfactory  for  stability  and 
control  calculations.  An  approach  such  as  developed  in 
Reference  Ik  is  too  complex  for  use  here.  Hie  approximation 
suggested  in  Reference  2  appears  to  be  incorrect  since  it 
assumes  that  there  is  no  change  in  wing  downvash  field  and 
only  a  change  in  the  mass  flow  influenced  by  the  wing. 
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However,  If  there  is  no  change  in  the  dovnwash  behind  the 
wing,  then  from  the  viewpoint  of  vortex  theory  there  would 
be  no  change  in  wing  lift,  contradicting  the  momentum  theory 
result.  It  is  important  to  note  that  momentum  theory  cannot 
be  used  alone  and  must  be  combined  with  a  sectional  approach 
to  obtain  a  defined  answer.  A  more  recent  analysis,  pre- 
sented  in  Reference  16,  also  seems  to  be  incorrect  nnd  there¬ 
fore  will  not  be  used.  In  this  reference,  the  boundary 
conditions  on  the  problem  are  not  properly  satisfied. 

•Rierefore,  since  precise  theoretical  results  are  not  particu¬ 
larly  satisfactory,  it  is  considered  that  a  semiempirical 
approach  to  the  problem  that  agrees  with  the  various  limiting 
cases  involved  is  the  best  alternative.  As  a  result,  it  is 
assumed  that  the  lift  curve  slope  of  the  wing  based  on  re¬ 
sultant  velocity  in  the  slipstream  may  be  expressed  in  terms 
of  the  lift  curve  slope  of  the  wing  in  a  uniform  flow  as 

o 

V  cos(  i  +  a  -  a  ) 
w  e 


K  is  an  empirical  factor  that  may  be  adjusted  to  match  the 
low-speed  wing  characteristics  when  experimental  data  are 
available.  It  Is  assumed  that  the  variation  in  lift  curve 
slope  depends  upon  the  ratio  of  the  component  of  the  free- 
stream  velocity  parallel  to  the  slipstream  velocity 
[V  cos(iw  +  a  -  ae)]*  In  hovering  then,  when  V  approaches 

zero,  the  wing  lift  curve  slope  is 


(1  -  K) 


and  in  high-speed  forward  flight,  when  VR  =  V  cos(iw  +  a  -  ae)> 


The  empirical  constant  K  depends  in  a  complex  way  upon  the 
wing  geometry,  the  propeller  geometry,  and  the  immersed 
aspect  ratio,  as  well  as  other  factors.  To  obtain  good 
agreement  between  the  experimental  and  theoretical  values  of 
the  speed  stability  (f^)  in  hover,  a  value  of  K  =  0.5  was 
selected.  This  expression  will  be  discussed  further  in  the 
section  on  stability  derivatives. 

4.  The  angle  of  attack  of  the  wing  in  the  slipstream  ae  is 
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assumed  to  be  given  by  the  vector  sum  of  i'ully  developed  in¬ 
duced  velocity  produced  by  the  propeller,  2v,  and  the  free- 
stream  velocity,  V  (References  1,  2,  3  and  Figure  A-l). 

This  particular  aspect  of  the  wing  slipstream  problem  has  not 
been  considered  theoretically  in  detail  in  this  report  or  in 
other  analyses.  That  is,  all  of  the  above  analyses  have  in¬ 
vestigated,  precisely  speaking,  the  change  in  the  lift  of  a 
wing  immersed  in  a  slipstream  when  the  wing  geometric  angle 
of  attack  in  the  slipstream  is  changed,  rather  than  that  due 
to  a  change  in  the  direction  or  magnitude  of  the  external 
flow.  It  is  implied  that  these  two  sources  of  angle  of  attack 
changes  are  equivalent.  Therefore,  the  above  assumption  is 
made  in  absence  of  a  better  approach.  As  a  result,  the 
effective  wing  angle  of  attack  is  taken  to  be 


V  sin(i  +  a)  V  6in(i  +  a) 

-1  v  w  w  ,  , 

a  =  tan  x  — — •  — — —  ■■  — =  sin  ■  (ll) 

e  V  cos(i  +  cr)  +  2v  V_ 


w 


R 


and  the  following  expressions  are  used  for  wing  lift  and  drag. 


L  =  i  P  V2  S  C  (a  +  a  ) 

s  2  I!  L,S  o  e' 

e 


(12) 


CL,S 


Ds  -  I  0  4  S  [CD,S„  *  VaT  (*«  ’  “e'2]  <«> 


The  lift  is  taken  perpendicular  to  the  slipstream  velocity 
VR  and  the  drag  parallel  to  Vp  .  For  simplicity,  the  drag 
is  computed  on  the  basis  of  an  elliptically  loaded  wing,  while 
in  the  actual  case  the  loading  may  be  far  from  elliptical. 

Thus  the  propeller  induced  flow  has  been  Idealized  to  a  uniform  tubular 
flow  of  velocity  2v0  that  is  added  vectorially  to  the  free  stream.  The 
wing  is  assumed  to  be  completely  immersed  in  this  flow  and  only  the  finite 
size  of  the  slipstream  is  taken  into  account  in  the  computation  of  wing 
lift  and  drag.  These  assumptions  imply  that  the  design  and  relative  lo¬ 
cation  of  the  propeller  on  the  wing  have  no  influence  on  the  wing  lift. 

HORIZONTAL  TAIL 


The  horizontal  tail  lift  is  calculated  on  the  basis  of  a  simplified  approxi¬ 
mation  to  wing  dovnwash  assuming  an  elliptically  loaded  wing,  that  is. 
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It  is  assumed  that  the  downwash  velocity  is  the  same  inside  and  outside  the 
slipstream.  The  following  expression  is  then  obtained  for  the  downwash 
angle  outside  the  slipstream 


e  =  -2-  C  (l  _  j£)  HE  (q-  +  a  )  (15) 

tt  A  L,Sa  T;  v  o  e' 

The  downwash  reduction  factor,  M^,  accounts  for  the  location  of  the 

horizontal  tail  with  respect  to  the  wing  trailing  vertex  system.  If  the 
distance  of  the  horizontal  tail  above  the  wing  wake  is  small  compared  to 
the  semispan  of  the  wing,  then  Mp  may  be  expressed  as 


„  <HT 

T  b 

w 
2 

The  tail  efficiency  T)T  is  assumed  to  be  one  (t^  =  1.0)  because  of  a 

lack  of  experimental  data.  To  obtain  accurate  information  on  the  down- 
wash,  experiments  are  necessary,  and  Equations  (l4)  through  (l6)  are  used 
only  to  obtain  an  estimate  of  the  tail  contributions.  Note  that  experi¬ 
ments  should  include  a  measurement  of  the  velocity  at  the  tail  to  determine 
T),p,  as  well  as  measurement  of  the  downwash  angles. 


WING  FLAP  CHARACTERISTICS 


The  wing  flap  characteristics  used  are  shown  in  Figure  37.  The  changes  in 
C_  and  aQ  are  included  directly  in  Equations  (12)  and  ( 13 )  for  wing 
uo 

lift  and  drag.  The  pitching  moment  due  to  flap  deflection  is  assumed  to 
be  proportional  to  slipstream  velocity  squared,  and  independent  of  the 
ratio  of  slipstream  velocity  to  free-stream  velocity.  While  there  would 
be  an  apparent  camber  change  due  to  the  finite  extent  of  the  slipstream 
(Reference  15 ),  this  effect  is  assumed  to  be  small  in  the  absence  of  ex¬ 
perimental  data. 

FUSFLAGE  PITCHING  MOMENT  AND  DRAG 

Again,  for  lack  of  experimental  data,  the  fuselage  pitching  moment  charac¬ 
teristics  are  computed  on  the  basis  of  slender  body  theory  (Reference  11, 
page  64)  using  the  following  equation; 

Mf  =  2(Vol)f  Kf  q  a  (17) 
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The  fuselage  drag  is  expressed  in  terms  of  an  equivalent  flat  plate  drag 
area  f  (Reference  10,  page  221)  as 

Df  =  q  f  (18) 

This  equivalent  flat  plate  area  is  assumed  to  be  independent  of  angle  of 
attack  and  is  taken  to  be  5*0  square  feet.  The  frontal  area  of  the 
fuselage  is  92  square  feet  and  the  fineness  ratio  is  5*3*  Assuming  that 
the  fuselage  is  a  streamline  body,  from  Reference  21,  page  6-19,  Figure 
25,  a  drag  coefficient  based  on  frontal  area  of  0.054  is  obtained,  giving 
the  above  value  of  f.  Experimental  data  are  necessary  to  obtain  a  better 
estimate  of  the  fuselage  characteristics. 

Further  experimental  data  are  desirable  to  obtain  better  estimates  of  the 
two  parameters  and  K.  Since  the  analysis  indicates  that  the  speed 

stability  (Mu)  is  dominated  by  the  wing  contribution  in  hover,  it  is 

considered  that  the  empirical  factor  K  may  be  estimated  reasonably  well 
by  use  of  test  results  near  hover.  Wake  surveys  are  desirable  to  obtain 
substantiation  of  the  assumed  value  of  T^,  and  the  value  of  one  should 

be  considered  only  as  a  rough  estimate. 
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DISCUSSION  OF  RESULTS 


TRIM  CONDITIONS 


The  force  equations  for  level  flight  equilibrium,  parallel  and  perpen¬ 
dicular  to  the  propeller  shaft  (Figure  A-l)  are 


T  +  L„  sin  0/  +  D„  cos  a  +  D~  cos  i  =  W  sin  i  (19) 

s  e  3  e  f  w  w  '  ' 


L,.  cos  c*e  +  sin  +  Df  sin  iy  +  H  =  W  cos  i^  (20) 

It  is  assumed  that  forces  produced  by  the  horizontal  tail  and  the  tail 
rotor  do  not  influence  the  balance  of  forces  significantly,  so  that 
Equations  (l))  and  (°0)  may  be  solved  for  the  flight  condition.  The 
moment  equation  will  determine  the  longitudinal  control  setting  necessary 
for  pitching  moment  trim. 

The  quantities  on  the  left-hand  side  of  Equations  (19)  and  (20)  are 
functions  of  the  following  variables,  in  addition  to  the  airplane  geometry: 

0  -  propeller  blade  angle  (radians) 

0  •  ( 

V  -  flight  velocity  (feet  per  second) 

Q  -  propeller  rotational  speed  (radians  per  second) 

6j.  -  flap  deflection  (radians  or  degrees) 

v  -  propeller  induced  velocity  (feet  per  second) 
a  -  aircraft  angle  of  attack  (radians  or  degrees) 
i  -  wing  incidence  angle  (radians  or  degrees) 

The  momentum  equation  relates  propeller  thrust  to  induced  velocity.  As  a 
result,  there  are  three  equations  [(l),  (l?)  and  (20)]  with  seven  unknowns. 
The  wing  incidence,  flap  deflection  (flap  deflection  as  a  function  of  wing 
incidence  is  shown  in  Figure  A-2),  propeller  rotational  speed,  and  air¬ 
craft  angle  of  attack  (zero)  are  specified.  Then  the  propeller  blade 
angle,  trim  velocity,  and  propeller  induced  velocity  can  be  determined. 

The  results  of  this  calculation  are  shown  in  Figures  1  and  3  as  a  function 
of  wing  incidence  angle.  Figure  2  presents  propeller  thrust  and  wing  lift 
as  a  function  of  incidence  angle.  The  slipstream  velocity  and  the  ef¬ 
fective  wing  angle  of  attack  are  shown  in  Figures  4  and  5* 
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In  order  to  obtain  an  estimate  of  the  manner  in  which  various  assumptions 
enter  into  the  equilibrium  calculations,  a  first-order  approximation  to 
Equations  (19)  and  (20)  may  be  obtained  by  neglecting  all  left-hand  terms 
except  thrust  in  Equation  (19),  and  the  wing  lilt  in  Equation  (20): 


T  ^  W  sin  iw 


(19a) 


Ls  cos  Og  —  W  cos  iy 


(20a) 


It  can  be  seen  from  Figure  2  that  Equation  (19&)  gives  a  very  good  approxi 
mation  to  the  thrust  required  for  equilibrium. 

The  wing  lift  in  terms  of  flight  variables  is  given  by  Equation  (12)  and 
(ye  by  Equation  (11).  Assuming  that  effective  wing  angle  of  attack  is 

small,  Equation  (20a)  may  be  expressed  as 


P  2  r,  x  V  C0S(iw4«-eire)T 

Ls  .  |  *  S  CL  [(1  -  K)  ♦  K  - - - I»c 


V  Bin(i  ■*) 

- — 


W  cos  i 


w 


cos  (V 


(20b) 

Since  the  slipstream  velocity,  Vp,  is  roughly  constant  through  transition 

(Figure  4),  Equation  (20b)  may  be  viewed  as  primarily  determining  the  re¬ 
lationship  between  flight  velocity  and  wing  incidence  angle.  At  a  given 
wing  incidence,  changing  aQ  by  deflecting  a  flap  will  have  a  strong 

influence  on  the  trim  velocity  of  the  aircraft. 


The  influence  of  flap  effectiveness 


fop 

96f 


on  trim  velocity  is  shown  in 


Figure  1.  A  plus  or  minus  33  l/3-percent  change  in  flap  effectiveness 
causes  a  large  variation  in  the  relationship  between  wing  incidence  and 
flight  velocity  as  shown.  No  variation  occurs  at  wing  incidence  higher 
than  8cP  because  the  flap  is  not  deflected  (Figure  A-2). 


The  variation  of  slipstream  velocity  with  wing  incidence  is  comparatively 
small  because  as  the  flight  speed  increases,  the  induced  velocity  de¬ 
creases  (Figure  3),  resulting  in  a  slowly  varying  slipstream  velocity. 
Calculation  of  the  slipstream  velocity  is  not  particularly  sensiiive  to 
flap  effectiveness  (Figure  4). 

The  effective  wing  angle  of  attack,  ae,  particularly  the  maximum  value, 
is  sensitive  to  flap  effectiveness  (Figure  5)>  as  expected  from  the  above 
discussion.  Thus,  the  flap  may  be  viewed  as  a  means  of  controlling  the 
effective  wing  angle  of  attack  as  seen  from  Equation  (20b).  The  variation 
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of  the  parameters  in  the  expression  for  the  effective  wing  angle  of  attack 
is  such  that  a  maximum  value  of  ae  occurs  at  a  wing  incidence  of  approxi¬ 
mately  2^3  . 

The  pitching  moment  acting  on  the  airplane  with  no  tail  rotor  thrust  and 
no  horizontal  tail  lift  is  shown  in  Figure  6.  The  important  contributions 
to  the  pitching  moment  acting  on  the  airplane  from  the  various  components 
are  also  shown.  Ihe  moments  arising  from  wing  lift  and  propeller  thrust 
vary  with  wing  incidence  because  of  changes  in  the  magnitudes  of  the  forces 
and  the  moment  arms  about  which  they  act.  An  increasing  nose-up  pitching 
moment  from  the  sum  of  the  thrust  and  wing  lift  with  decreasing  wing  inci¬ 
dence  is  typical  of  tilt-wing  configurations  for  normal  center- of-gravity 
locations  (Reference  18).  The  contribution  of  the  propeller  pitching 
moment  due  to  the  induced  velocity  distribution  over  the  propeller  disc  is 
small.  The  flap  characteristics  have  a  significant  influence  on  the  moment 
required.  To  obtain  an  accurate  prediction  of  the  pitching  moment  required 
to  trim,  precise  information  on  the  pitching  moments  produced  by  flap  de¬ 
flection  in  a  slipstream  is  necessary. 

In  conclusion,  it  may  be  noted  that  the  magnitudes  of  the  propeller  thrust 
and  wing  lift  are  primarily  determined  by  wing  incidence  angle.  The  trim 
velocity  corresponding  to  this  incidence  angle  depends  strongly  on  the 
assumed  relationship  between  the  wing  lift  and  the  flight  variables. 
Computations  of  both  trim  velocity  and  pitching  moment  are  sensitive  to 
flap  characteristics.  To  obtain  a  good  prediction  of  the  wing  incidence/ 
trim  velocity  relationship,  it  is  necessary  to  have  good  estimates  of  the 
wing  lift  curve  slope  and  flap  effectiveness  of  a  wing  immersed  in  a  high 
velocity  slipstream. 

STABILITY  DERIVATIVES 


The  stability  derivatives  of  the  example  airplane, with  respect  to  a  sta¬ 
bility  axis  system  for  various  level  flight  equilibrium  conditions,  were 
calculated  using  the  theoretical  models  discussed  in  the  first  section. 
This  section  presents  a  detailed  discussion  of  the  trends  of  the  deriva¬ 
tives  and  the  contributions  of  the  airplane  components  to  each  of  the 
derivatives.  Detailed  derivations  of  various  formulae  discussed  are  given 
in  Appendix  B. 

Force  Derivatives  With  Respect  to  Horizontal  Velocity  (\  ,  Z^) 


The  force  derivatives  arise  primarily  from  propeller  thrust  and  wing  lift 
variations.  Recall  that  wing  lift  and  drag  are  taken  perpendicular  and 
parallel  to  the  direction  of  slipstream  velocity.  Generally,  it  was  found 
that  the  contributions  arising  from  the  propeller  inplane  force  and  the 
wing  drag  were  small.  Therefore,  only  the  rate  of  change  of  thrust  and 

wing  lift  with  horizontal  velocity  (  ^  will  be  discussed  in 

V  du  5u  / 

relation  to  the  airplane  derivatives. 

It  is  necessary  first  to  consider  the  rate  of  change  of  effective  wing 
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angle 


du 


of  attack  and  slipstream  velocity  with  horizontal  velocity 
dVR  N 

,  —  ]  as  well  as  the  changes  in  propeller  induced  velocity, 
du  J 


Analytically,  it  is  more  convenient  to  consider  the  propeller  induced  ve¬ 
locity  derivatives  with  respect  to  the  nondinensional  velocities  u  and 
p,  ,  perpendicular  and  parallel,  respectively,  to  the  propeller  plane.  In 
fact,  it  was  noted  during  the  course  of  this  analysis  that  the  aerodynamic 
force  and  moment  expressions  became  considerably  simplified  when  treated 
in  terms  of  the  variables  u  and  p.  with  an  axis  system  aligned  with  the 
propeller  shaft.  This  approach  has  the  disadvantage  that  the  inertia  and 
gravity  terms  will  appear  in  an  unconventional  form  in  the  equations  of 
motion,  so  it  was  decided  to  proceed  in  a  more  conventional  manner  and 
accept  the  resulting  complication  in  the  aerodynamic  terms. 

First,  the  rate  of  change  of  propeller  induced  velocity  with  a  change  in 

*Xi 

velocity  parallel  to  the  propeller  plane  ±  is  considered.  The  ex- 

dn 

press  ion  as  given  in  Appendix  B  is 


where 


*  (2D 


2 

+  » 


The  first  term  in  the  numerator  of  Equation  (21)  arises  from  the  thrust 
variation  due  to  the  average  increase  in  the  free-stream  velocity  component 
over  the  disc,  and  the  second  term  accounts  for  the  variation  in  induced 
velocity  resulting  from  the  change  in  mass  flow.  The  second  term  generally 
dominates  this  expression  and  is  roughly  proportional  to  X^  p,  ,  since  § 
is  almost  constant  through  transition  as  shown  in  Figure  8.  Since  p.  is 
VR 

equal  to  -£  sin  a„  ,  p,  varies  approximately  as  cr_  .  At  high  wing  inci- 
QR  e 

dences,  pi  increases  more  rapidly  than  the  induced  velocity,  X^ 
decreases  (with  a  change  in  flight  velocity),  and  therefore  the  derivative 
3X 

— ±  increases  in  value.  At  low  wing  incidences,  both  X  and  p  de- 

dp.  1 

crease  with  increasing  flight  velocity,  and  so  this  derivative  decreases. 

The  variation  of  this  derivative  with  wing  incidence  will  be  similar  to 
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-*Y 


%  ^ 


A 

\ 


the  variation  of  a 

ax 


Although  there  is  some  influence  on  the  derivative 


W  A  • 

— 1  from  the  denominator,  the  major  variation  arises  from  the  numerator. 

ax4 


The  trends  of  the  derivative 


‘i 


b\i 


are  shown  in  Figure  7*  Physically,  the 


magnitude  of  this  derivative  indicates  that  for  the  velocity-propeller 
angle  of  attack  relationship  encountered  during  a  transition,  the  propeller 
induced  velocity  is  rather  insensitive  to  velocity  perturbations  parallel 
to  the  plane  of  the  propeller. 

ax. 

Trie  other  induced  velocity  derivative,  —1  ,  may  be  expressed  as  ■ 

av 


1  racr 

1  L8 


1  + 


1  T  an 

s  Lr 


X  ^  +  v 

r 

X^  +  v 


]] 

)] 


(22) 


Again,  since  £  is  roughly  constant  through  transition,  the  variation  of 
this  quantity  is  dominated  by  the  change  in  Lnduced  velocity,  , 

X^  +  v 

through  the  transition,  since  -  ■  ■  =  1.0  (Figure  8).  Also  because  the 

? 

resultant  velocity  at  the  propeller  plane,  5  >  Is  approximately  constant 
throughout  most  of  the  transition  (Figure  8),  the  induced  velocity  varies 
as  the  thrust  [see  Equation  (l)].  This  derivative  is  shown  as  a  function 
of  wing  incidence  angle  in  Figure  7-  Physically,  the  magnitude  of  this 
term  indicates  that  variation  in  propeller  induced  velocity  acts  to  reduce 
greatly  the  effect  of  a  velocity  perturbation  normal  to  the  propeller 
plane  on  the  slipstream  velocity. 

The  variation  in  the  slipstream  velocity  (see  Appendix  B)  is  given  by 


ax 

i  2\.  +  v  s  SX.  2(  2X .  +  v)  y—  +  pi 

—  — -  =  — i - (  1  +  2  — £  )  cos(i  +  cy)  +  - - - sin(i  +  a) 

qr  au  Vp  v  av  y  v  w  vn  v  w  ' 


(23) 


The  first  term  represents  the  effect  Ox  the  component  of  horizontal  ve¬ 
locity  normal  to  the  propeller  plane,  and  the  second  term  is  due  to  the 
component  of  horizontal  velocity  parallel  to  the  propeller  plane. 


Since  the  direction  of  the  resultant  slipstream  velocity  is  always  closely 
aligned  to  the  propeller  shaft  direction  (the  effective  angle  of  attack  is 
small),  a  change  in  the  free-stream  velocity  component  perpendicular  to  the 
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propeller  plane  has  only  a  small  influence  on  the  magnitude  of  the  re¬ 
sultant  velocity.  As  discussed  earlier,  the  induced  velocity  change 
arising  from  thi;  component  of  free-stream  velocity  is  small.  The  first 
term  in  Equation  (23),  arising  from  the  change  in  the  free-stream  velocity 
component  along  the  propeller  shaft,  will  dominate  this  expression.  How¬ 
ever,  this  term  will  also  be  small  due  to  the  compensating  effect  of  pro¬ 
peller  induced  velocity  variation  with  v  discussed  above.  In  hovering, 


(^1  +  2  — i  ^  is  a  small  negative  quantity.  As  the  flight  velocity  in¬ 
creases,  I — H  decreases  and  fl  +  2  — -  ^  becomes  zero  and  then  a  small 

ov  V  av  J  9V 

positive  quantity  with  further  increase  in  trim  velocity.  Therefore,  — - 

s . 3VR  ,  5u 

varies  as  shown  in  Figure  10,  and  is  a  small  quantity  (^|--£|  <  0.2  J  . 


The  result  that  this  term  is  small  will  usually  be  the  case  as  seen  from 
the  following.  In  hovering,  |  =  X^  ,  and  therefore,  from  Equation  (22) 


(22a) 


ax. 

Thus,  in  hovering,  -  0.5  <  — =  <  -  1.0  .  As  the  velocity  component  normal 

dv 

to  the  propeller  plane  increases,  the  induced  velocity,  X^  ,  will  de¬ 
crease,  and  also  its  rate  of  change  with  v  will  decrease. 


The  magnitude  of  the  term  ^1  +  2  g— i  ^  is  very  sensitive  to  the  as¬ 
sumption  that  the  slipstream  velocity  is  computed  on  the  basis  of  the 
fully  developed  induced  velocity.  It  is  possible  that  a  factor  less  than 

ax. 

two  should  appear  in  front  of  —  to  account  for  the  actual  variation  in 

3v 

induced  velocity  across  the  wing  chord. 


As  the  trim  velocity  corresponding  to  a  given  wing  incidence  becomes 

avR 

smaller,  the  larger  the  minimum  value  of  — -  will  be.  Figure  9  shows 

3u 


avR 

lines  of  constant  — -  as  a  function  of  trim  speed  and  wing  incidence 

3u 

angle.  The  upper  boundary,  for  which  — s  equal  to  0.2,  occurs  at  too 

3u 

large  a  velocity  to  be  shown  on  this  graph.  It  therefore  appears  that  in 
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the  typical  situation,  the  slipstream  velocity  will  be  relatively  insensi¬ 
tive  to  changes  in  horizontal  velocity  if  the  assumption  regarding  the  full 
development  of  the  slipstream  is  valid. 

The  fourth  quantity  of  interest  here,  the  rate  of  change  of  wing  effective 
angle  of  attack  with  horizontal  velocity  (see  Appendix  B)  may  be  expressed 
as 


^e  _  sln(iw  +  a)  A  V.  \ 

au  vR  V.  yR  au  / 


Since,  in  this  example,  the  trim  velocity  always  lies  well  within  the 

I  <  0.2^  and  i  is  small  when  the  trim  velocity  is  within 
J  VR 

this  boundary,  a  good  approximation  to  Equation  (24)  is  obtained  by  ne¬ 
glecting  the  second  term 


avR 

boundary  (  — C 

vau 


aoe  _  sin(iw  +  Of) 
”  =  VR 


(25) 


For  the  example  configuration  where  a  is  equal  to  zero,  the  agreement 
between  Equations  (24)  and  (25)  is  shown  in  Figure  11.  Physically,  this 
approximation  implies  that  the  primary  change  in  effective  wing  angle  of 
attack  due  to  a  horizontal  velocity  perturbation  arises  from  the  change 
in  the  component  of  the  horizontal  velocity  parallel  to  the  propeller 
plane,  and  the  effects  of  component  normal  to  the  propeller  plane  and  the 
variation  in  induced  velocity  are  small.  As  the  wing  incidence  is  reduced, 
the  increasing  effect  of  the  component  parallel  to  the  shaft  is  compensated 
by  the  induced  velocity  change. 

The  rate  of  change  of  propeller  thrust  with  horizontal  velocity  may  be 
expressed  as 


ai 

}U 


=  p  H  R2  (QR)  92  [-  i  (l  +  !ii  )  cos(iw  +  a)  + 


C3e0.75R"-2g“)sin(iw+^] 


(26) 


16 


.  — ■ -X 


v? 


i- 


This  derivative  is  shown  as  a  function  of  incidence  angle  in  Figure  12. 
The  first  term  dominates  and  causes  the  value  to  increase  (approximately) 
as  cos( 1  +  a) . 

The  rate  of  change  of  wing  lift  with  horizontal  velocity  may  be  expressed 
as 


du 


±  p  S  C  „ 
2  H  L,S 


a , 


2 

VZ  — -  +  p  s  c  „ 

R  au  H  l,s 


(°o  +  ae') 


+ 


1 

2 


S  V, 


ac 


+  ae^ 


L,S 


a. 


&u 


(27) 

Thus,  variations  in  wing  lift  with  horizontal  velocity  may  be  considered 
as  arising  from  three  sources:  first,  the  change  in  effective  angle  of 
attack;  second,  the  change  in  slipstream  velocity;  and  third,  the  change 
: n  wing  lift  curve  slope.  The  variation  in  this  derivative  with  incidence 
angle  and  the  contribution  of  each  of  the  above  terms  are  shown  in  Figure 
13.  The  first  term  dominates  at  high  incidence  angles,  and  the  last  terra 
becomes  quite  important  at  low  incidence  angles.  The  last  term  must  be 
considered  questionable  due  to  lack  of  knowledge  of  the  exact  functional 
relationship  for  the  lift  curve  slope.  This  particular  effect  was  not 
taken  into  account  in  either  Reference  3  or  its  equivalent  form  in  Refer¬ 
ence  2.  The  middle  term  seems  unimportant;  however,  it  would  be  larger  if 
the  slipstream  velocity  were  computed  on  the  basis  of  other  than  the  fully 
developed  induced  velocity.  As  a  result  of  these  various  contributions, 
this  derivative  remains  roughly  constant  throughout  the  transition. 

The  total  value  (summation  of  all  contributions)  of  these  two  stability 
derivatives  and  the  manner  in  which  they  vary  are  discussed  next.  First, 
consider  the  rate  of  change  of  horizontal  force  with  horizontal  velocity: 

♦  ran 

du  fuselage 


)  +  r  an 

v  du  ^fuselage 


dx  /  ax\  +r-^) 

du  V  du  ^propeller  ^  du  Aring 


or 


dL 


CLX  =  (  A  cos  i  -  (  — £  A  sin(i  -  < 

iu  Uu  J  w  v  au  J  v  * 


The  wing  and  propeller  contributions  to  this  derivative  are  shown  in 
Figure  14.  The  wing  contribution  dominates  at  high  incidenc^&nd  pro¬ 
peller  contribution  dominates  at  low  incidence.  At  high  inoapnce  angles. 
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s  e  r  '  w  ~\ 

— —  arises  primarily  from  —  =  |  — — — — —  as  shown  In  Figure  13. 

■  L  Vr  J 


au 


au 


The  rate  of  change  of  vertical  force  with  horizontal  velocity  may  be  ex¬ 
pressed  as 


az  =  (  12  n 

Su  au  / 


(¥-) 

v  au  wing 


au  /rotor 


or 


12 

au 


.  SL  a® 

^  sin  itt  -  —2  cos(iw  -  aj  -  LR  sin(  i„  -  arp)  ^—2 


au 


v  9u 


w 


du 


(29) 


The  behavior  of  this  derivative  with  wing  incidence  is  shown  in  Figure  15. 

The  derivative  is  dominated  by  the  wing  lift  contribution  — -  cos(i  -  a  ) 

au  w 

throughout  the  transition  regime. 


Force  Derivatives  With  Respect  to  Vertical  Velocity  ( X v  ,  Z^.) 

First,  the  variations  in  effective  wing  angle  of  attack,  slipstream  ve¬ 
locity,  propeller  thrust,  ano  wing  lift  with  vertical  velocity  are  dis¬ 
cussed,  and  then  the  complete  vehicle  derivatives. 

First,  the  slipstream  velocity  derivative  may  be  expressed  as 


av_  2X  +  u  dX 

A  -2  -f  — - V  1  .  2  _ly  .  sin(l  ♦*)> 

QR  aw  V  VR  /  V  Sy  /  V  w  / 


ax. 

2(2X  +  V)  -i 

i  dpt 


+  \i 


cos( i  +  a) 


(30) 


The  variation  of  this  quantity  is  shown  in  Figure  l6.  At  high  incidence 

ax. 

angles  the  first  term  dominates.  Since  both  sin(  i  +  cr)  and  (l  +  ?  — — ) 

w  ou 

decrease  with  decreasing  wing  incidence,  the  first  term  decreases  and  at 
wing  Incidence  angles  below  30°  the  second  term  becomes  important.  The 
rate  of  change  of  effective  wing  angle  of  attack  with  vertical  velocity  is 
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dv 


(31) 


C0S(1W  +  O')  (1 
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Of 


cos( 


fp 

dw 


dw  \ 

iW  +  ®) 


For  the  example  configuration,  where  the  fuselage  angle  of  attack  is  equal 

Sop 

to  zero,  Figure  17  indicates  that  the  quantity  — —  is  dominated  by  the 

dw 

first  term  and,  similarly  to  — =,  arises  primarily  from  the  change  in  ve- 

du 

locity  component  nonnal  to  the  wing;  it  is  not  significantly  influenced  by 
the  second  term.  The  rotor  induced  velocity  variation  is  compensated  for 
by  the  velocity  component  variation  parallel  to  the  shaft.  Again  the 
magnitude  of  this  compensating  effect  depends  upon  the  assumption  of  a 
fully  developed  induced  velocity. 

The  second  term  of  Equation  (31)  is  a  rather  strong  function  of  the  trim 
velocity  due  to  the  appearance  of  the  effective  wing  angle  of  attack.  In 
Figure  9>  showing  trim  velocity  plotted  versus  wing  incidence  angle,  the 
boundaries  where  the  second  term  is  10  percent  and  20  percent  of  the  first 
term  are  shown.  As  can  De  seen,  it  is  possible  that  at  high  incidence 
angles,  the  second  term  may  become  significant. 

The  rate  of  change  of  thrust  with  vertical  velocity  can  be  expressed  as 
-  =  p  n  EL  (QR)  —  \(l  +  —  ^  sin(iw  +  q0  + 

dw  2  LV  3v  J 

C-  -^)C0E(lw  +  a)] 

(32) 

This  quantity,  shown  in  Figure  18,  again  at  high  incidence  angle  is 
dominated  by  the  first  term.  Ihe  second  term  increases  with  decreasing 
incidence,  becoming  important  at  low  tilt  angles.  As  a  result,  this  de¬ 
rivative  remains  roughly  constant  about  equal  to  its  hovering  value.  In 
hovering, 
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From  Equation  (22),  with  ?  =  and  v  =  0, 


axi  C  f  +  vi  ) 


dv 


(f  +  M 


and  from  Equation  (l),  with  v  =  g.  =  0, 
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Thus,  since  —  will  be  relatively  independent  of  the  configuration  of  the 
a 

vehicle  (Reference  10,  page  64),  this  derivative  will  primarily  vary  with 

Crp 

tip  speed,  increasing  as  the  tip  speed  decreases  at  constant  —  . 

The  wing  lift  variation  with  vertical  velocity  is  given  by 
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(35) 


Similarly,  three  contributions  appear  as  discussed  with  respect  to  Equation 

dL 

(27).  The  numerical  values  of  — -  are  plotted  in  Figure  19 .  The  first 

&W  d<*e 

term  dominates  throughout  transition.  The  term  —  is  nearly  pro- 

d* 

portioned,  to  c^s(iw  +  a)  as  noted  in  the  discussion  of  Equation  (31). 
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Physically,  as  the  wing  incidence  decreases,  the  variation  in  the  component 
of  vertical  velocity  normal  to  the  propeller  shaft  increases,  causing  a 
larger  effective  angle  of  attack  change  from  the  same  change  in  vertical 
velocity.  Again,  the  last  term  should  be  considered  rather  questionable 
due  to  a  lack  of  precise  knowledge  of  the  wing  lift  curve  slope. 

The  total  value  of  the  airplane  derivatives  (summation  of  all  contri¬ 
butions)  is  now  considered.  First,  the  rate  of  change  of  horizontal  force 
with  vertical  velocity  is 
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This  derivative  is  shown  in  Figure  20.  At  high  tilt  angles,  £ 

aw 

aT  *LS 

sin(L.  -  a  )  and  2_  cos  i  dominate.  Since  - —  sin(i„  -  ff„)  1b 

W  e  gW  W  QW  w  t. 

aT  ax 

greater  than  ~  cos  i„  at  high  tilt  angles,  2-  is  negative.  As  the 
aw  aw 

/  X 

incidence  decreases,  Ls  cos(iw  -  ae)  — -  becomes  increasingly  important 


3w 


aL„ 


■*Y  **  g 

and  —  becomes  positive.  The  term  involving  —  is  due  to  the  change 

aw  aw 

aofe 

in  magnitude  of  the  wing  lift,  and  the  term  Ls  cos(iw  -  ae)  —  is  due  to 

aw 

the  tilt  of  the  wing  lift  vector.  At  high  incidence  angles  the  change  in 
lift  is  important  and  at  low  incidence  angles  the  tilt  of  the  lift  vector 
becomes  significant,  since  the  magnitude  of  the  wing  lift  increases  as  the 
wing  incidence  deci'eases  (Figure  2). 


It  was  noted  from  experimental  measurements  that  the  variation  of  the 
horizontal  force  with  vertical  velocity  is  quite  nonlinear.  Bie  theo¬ 
retical  computation  of  the  horizontal  force  as  a  function  of  angle  of 
attack  using  the  expression  given  in  Appendix  B  is  shown  in  Figure  38. 

The  nonlinear  variation  in  this  force  arises  from  the  change  in  the  magni¬ 
tude  of  the  lift  vector  as  it  tilts.  That  is,  as  the  above  term, 

do 

L  cos(i  -  a  )  — ,  becomes  important,  the  wing  lift  also  varies 
s  w  e  Bw 
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appreciably  with  effective  wing  angle  of  attack,  resulting  in  this  tern's 
increasing  with  <ye  .  This  is  probably  not  a  particularly  important 

effect  on  the  dynamics,  since  this  derivative  itself  usually  is  not  im¬ 
portant.  The  nonlinearity  is  particularly  noticeable  where  the  linearized 
derivative  is  near  zero  as  shown  in  Figure  38. 

The  rate  of  change  of  vertical  force  with  vertical  velocity  may  be  ex¬ 
pressed  as 


or 


3Z  _  aZ 

dw  dw  -  propeller 


sin(iw  -  ae) 


(37) 


The  behavior  of  this  derivative  is  shown 


in  Figure  21. 


At  large  incidence 


angles , 


~  sin  i„  dominates.  As  the  wing  incidence  decreases,  —  sin  i,, 
9w  Sw 


3Ls 

decreases  and  - —  cos(iu  -  a„) 
aw  '  w  e' 

creases  roughly  proportional  to 

bution. 


becomes  dominant.  This  derivative  in- 
cos(iw  -  <*e)  due  to  the  wing  lift  contri- 


Pitching  Moment  Derivatives  With  Respect  to  Horizontal  Velocity,  Vertical 
Velocity  and  Pitch  Pate  ,  My  ,  ) 


The  pitching  moment  derivatives  are  now  considered.  First,  examine  the 
variation  of  the  pitching  moment  with  horizontal  velocity,  which  may  be 
broken  'down  into  the  following  contributions : 


aM  , 

'  aM  - 
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tail 

rotor 

(38) 
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As  may  be  seen  from  Figure  22,  the  primary  contributions  to  —  arise  from 

au 

the  propeller  and  the  wing.  The  major  contribution  from  the  wing  arises 
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from  the  term 


V%  S'  S  \ 

—  f  — —  cos  <y„  )  as  shown.  Since  l  decreases  roughly 

I  Uu  e  J  W 

in  proportion  to  sin  iw  (Figure  A-3),  the  wing  contribution  varies 

roughly  as  sin  iw  .  The  propeller  contribution,  arising  from  the  vari¬ 
ation  in  induced  velocity  across  the  propeller  disc,  Equation  (2),  has  a 
noticeable  contribution  in  hovering.  As  the  wing  incidence  decreases,  the 
propeller  contribution  decreases.  As  a  result  M,^  decreases  rapidly  as 

wing  incidence  is  decreased. 


The  contribution  of  the  horizontal  tail  to  this  derivative  will  depend 
upon  the  tail  incidence  and  downwash  characteristics  at  the  horizontal 
tail.  The  tail  contribution  may  be  evaluated  from  Equation  (l4). 
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T  HT 
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r  *  +  1^ 

du  du  J 


(39) 


The  first  term  depends  upon  the  tail  incidence  and  the  downwash  angle.  To 
estimate  the  importance  oi  this  term,  its  magnitude  for  a  tail  lift  coef¬ 
ficient  equal  to  1.0  is  shown  in  Figure  23.  While  the  actual  tail  lift 
coefficient  would  vary  in  some  complex  way,  this  indicates  the  importance 
of  the  term.  Thus,  as  a  result  of  the  tail  incidence,  there  may  be  a  sig¬ 
nificant  contribution  to  the  variation  of  pitching  moment  with  horizontal 
speed  at  lower  wing  incidence  angles.  An  increase  in  tail  efficiency, 

7]^,  would,  of  course,  increase  the  magnitude  of  this  term. 


The  second  term  arises  from  the  wing  and  propeller  downwash  variation  with 
horizontal  velocity.  The  downwash  behind  a  wing  immersed  in  a  high  ve¬ 
locity  slipstream  is  difficult  to  estimate,  and  so  only  the  following  rough 
indications  of  the  size  of  the  second  term  in  Equation  (39)  are  considered. 
The  propeller  downwash  contribution  would  be 


du  9u 


and  the  wing  downwash  would  be 


7u%-7bc^  C 1-^)y  7s 


m 


(41) 


These  two  terms  are  plotted  in  Figure  23.  Thus,  assuming  that  the  hori¬ 
zontal  tail  characteristics  are  proportional  to  free-stream  dynamic 
pressure,  rather  small  contributions  to  the  derivative  are  indicated. 

Further  experimental  data  are  necessary  before  downwash  effects  can  be 
evaluated  in  more  detail. 
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The  rate  of  change  of  pitching  moment  with  vertical  velocity  (attitude  sta- 


dM 


shown 

in  Figure 

24a, 

is  made  up  of  the  following  contributions 

)  +( 

dM  n, 

>  + 

V  dw 

^wing  ^ 

dw  / 

propeller  ^  ^horizontal 

(IT; 


tail 

rotor 


dw  / 


The  pitching  moment  derivative  contribution  of  the  fuselage, 

,  was  not  calculated  due  to  lack  of  experimental  data  and 

/fuselage 

is  not  included  in  this  summation. 

At  high  incidence  angles,  the  wing  and  propeller  terms  dominate.  A  break¬ 
down  of  the  various  wing  contributions  is  also  shown  in  Figure  24b.  As 
the  wing  incidence  decreases  and  the  trim  velocity  increases,  the  hori¬ 
zontal  tail  provides  the  dominant  tern,  as  in  a  conventional  aircraft.  At 
high  incidence,  that  is,  low  speed,  the  wing  and  rotor  combine  to  produce 
a  statically  unstable  aircraft  (My  positive).  As  the  flight  speed  is  in¬ 
creased,  the  effect  of  the  horizontal  tail  increases  and  the  aircraft  be¬ 


comes  statically  stable.  The  wing  incidence  angle  at  which 


aw 


has  a 


maximum  positive  value  is  strongly  dependent  upon  the  trim  velocity  corre¬ 
sponding  to  this  incidence,  due  to  the  increasing  importance  of  the  hori¬ 
zontal  tail.  The  horizontal  tail  contributions  have  been  computed  es¬ 
sentially  neglecting  the  effects  of  slipstream  velocity.  While  this  may 
affect  the  precise  numerical  values  obtained,  it  would  be  expected  that  the 
trend  of  Mw  would  not  be  significantly  altered.  This  theoretical  vari¬ 
ation  of  angle  of  attack  stability  agrees  with  the  general  behavior  of  this 
derivative  noted  in  Reference  9  for  various  VTOL  aircraft. 

During  the  course  of  these  computations,  it  was  also  noted  that  there  were 
nonlinearities  in  this  derivative.  Figure  39  shows  the  pitching  moment 
variation  versus  airplane  angle  of  attack  for  various  wing  incidences.  The 
nonlinearity  arises  primarily  from  the  change  in  direction  of  the  wing  lift 
with  vertical  velocity  and  is  particularly  noticeable  near  the  incidence 
angle  where  the  linearized  derivative  is  zero.  In  the  actual  case,  this 
derivative  may  be  considerably  more  nonlinear  than  indicated,  due  to 
fuselage  and  downwash  characteristics. 

Figure  25  shows  the  variation  of  the  pitch  damping  derivative  with  wing 
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Incidence.  Contributions  to  the  pitch  damping  of  the  vehicle  arises  from 
the  following  sources; 


+ 


propeller 


c  ' 

)  +( 

'  9M- 

) 

v  se . 

norizontal 

.  ae 

tail 

tail 

rotor 

At  high  incidence,  the  contributions  of  the  propeller  and  tail  rotor 
dominate  and  are  generally  rather  small.  As  the  forward  speed  increases, 
the  horizontal  tail  provides  an  increasing  damping.  The  wing  provides 
only  a  very  small  contribution,  as  would  be  expected.  It  should  be  re¬ 
called  that  the  propeller  damping  arises  primarily  from  the  change  in  in¬ 
flow  distribution  and  may  be  somewhat  overestimated,  since  the  change  in 
induced  velocity  has  not  been  included.  It  appears,  then,  that  any  sig¬ 
nificant  pitch  damping  arises  from  the  horizontal  tail,  and  that  the  rest 
of  the  components  give  generally  negligible  contributions. 

LOCUS  OF  ROOTS  MID  SENSITIVITY  TO  STABILITY  DERIVATIVE  VARIATIONS 


Using  the  values  of  the  stability  derivatives  calculated  in  the  previous 
section,  the  characteristic  roots  of  the  aircraft  have  been  calculated 
[(using  Equations  (B-4),  (B-5),  and  (B-6)]  and  are  shown  in  Figure  26  as 
a  function  of  wing  incidence.  Hie  horizontal  tail  contribution  to  My 
has  been  neglected. 

At  hovering  and  slew  forward  speeds,  the  dynamics  are  similar  to  a  con¬ 
ventional  helicopter  (Reference  ll).  The  modes  of  motion  consist  of  an  un¬ 
stable  oscillation,  a  real  convergence  that  is  relatively  fast,  and  a  slow 
convergence.  In  hovering,  the  vertical  degree  of  freedom  (vertical  ve¬ 
locity)  is  uncoupled  and  is  a  slow  convergence.  The  other  three  roots  com¬ 
prise  a  convergence  and  an  unstable  oscillation  due  to  the  coupling  of 
pitch  attitude  and  horizontal  velocity.  As  the  horizontal  velocity  in¬ 
creases  (wing  incidence  decreases),  the  period  of  unstable  oscillation  in¬ 
creases  and  the  oscillation  becomes  stable,  tending  towards  the  charac¬ 
teristics  of  a  conventional  airplane  phugoid  (Reference  11).  The  vertical 
motion  becomes  coupled  with  the  attitude  and  horizontal  velocity,  and  the 
two  real  roots  coalesce  and  then  break  away  from  the  negative  real  axis  to 
form  a  well  damped  motion,  similar  to  the  conventional  aircraft  short- 
period  motion  at  low  incidence  angles. 

Since  numerous  assumptions  have  been  made  in  the  foregoing  analysis,  it  is 
of  interest  to  evaluate  the  sensitivity  of  the  root  location  to  variations 
in  the  stability  derivatives.  In  order  to  evaluate  the  variations  in 
roots,  each  derivative  was  individually  reduced  to  zero  and  doubled.  Care 
should  be  taken  in  interpreting  these  loci  for  zero  values  of  the  important 
stability  derivatives,  since  nonlinearities  are  present  which  will  become 
more  important  when  the  linearized  derivative  is  zero.  Hie  resulting 
changes  in  the  characteristic  roots  of  the  vehicle  are  shown  in  Figures  27 
through  31. 
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tor  clarity,  the  roots  located  on  the  negative  real  axis  and  computed  for 
the  incremental  values  compared  to  the  basic  values  of  Figure  26  are  shown 
on  two  additional  "displaced"  real  axes  so  that  they  are  not  confusingly 
superposed,  thereby  making  it  impossible  to  discern  their  value  and  lo¬ 
cation.  Of  course,  they  are  located  on  the  abscissa  of  the  original  axis 
.ystem  at  the  appropriate  values  shown  on  the  displaced  axes.  Also,  in 
order  to  distinguish  between  the  two  oscillatory  modes  that  occur,  one 
mode  of  motion  is  referred  to  as  the  unstable  mode,  even  though  it  becomes 
table  at  higher  trim  velocities,  and  the  other  is  referred  to  as  the 
table  mode. 

Variation  in  Xu  (Figure  27) 


Xu  is  most  noticeable  in  influencing  the  damping  of  the  unstable  mode.  An 
increase  in  |XU|  increases  the  damping.  It  also  has  some  effect  on  the 

comparatively  rapid  convergence.  It  has  only  a  small  effect  on  the  stable 
mode  that  appears  at  lower  tilt  angles. 

Variation  in  Zy  (Figure  28) 


Increasing  the  magnitude  of  lift  curve  slope  of  the  airplane  Zy  increases 

the  damping  of  the  stable  mode  and  has  little  effect  on  the  unstable  mode. 

A  rather  radical  alteration  in  the  general  appearance  of  the  locus  occurs 
at  low  incidence  angles  when  Zy  is  set  equal  to  zero.  The  association  of 

the  unstable  mode  with  the  classical  phugoid  is  less  clear,  since  as  wing 
incidence  decreases,  this  mode  tends  towards  what  would  probably  be  con¬ 
sidered  the  short-period  motion.  Physically,  the  zero  lift  curve  slope 
would  represent  a  rather  unusual  airplane. 

Variation  in  (Figure  29) 


Increasing  M.  in  a  negative  sense  (stable)  at  high  incidence  angles  re- 
0 

duces  the  instability  of  the  unstable  mode  as  well  as  increases  the  damping 
of  the  stable  mode,  that  is,  it  contributes  favorably  to  both  modes  of 
motion. 

Variation  in  My  (Figure  30 ) 

At  high  incidence  angles,  an  increase  in  the  positive  value  (unstable)  of 
My  causes  a  loss  in  damping  of  the  unstable  mode.  At  low  incidence 

angles,  the  increasing  negative  value  (stable)  results  in  a  marked  increase 
in  the  damping  of  the  unstable  mode.  Increasing  the  stable  value  of  My 

at  low  incidence  angles  causes  a  marked  increase  in  the  frequency  of  the 
stable  mode,  as  would  be  expected. 


| 
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Variation  in  My  (Figure  31) 


Decreasing  My  causes  a  lengthening  of  the  period  and  a  slight  increase 
in  damping  of  the  unstable  mode.  Increasing  My  causes  a  decrease  in  the 
period  of  the  unstable  mode.  The  term  My  was  not  reduced  to  zero  since 
it  causes  a  radical  reorientation  of  tne  locus.  With  jMy  equal  to  zero 
in  hovering,  the  attitude  and  horizontal  velocity  are  essentially  un¬ 
coupled.  A  positive  value  of  this  derivative  is  typical  of  all  of  these 
vehicles  near  hovering. 

Variation  in  Zy  and  Xw 


The  variations  in  Zy  and  Xy  over  the  range  considered  were  found  to 
produce  only  small  changes  in  the  loci. 

GENERAL  TRENDS  OF  LOCI  OF  ROOTS 


In  order  to  understand  more  clearly  the  general  tendencies  of  these  loci, 
some  sketches  are  given  below  considering  only  the  effects  of  Mw  and 
My,  at  fixed  values  of  the  other  derivatives.  The  variations  of  these 
two  derivatives  dominate  the  general  shapes  of  the  loci  presented. 

In  the  case  where  M^,  =  0  and  My  is  varied,  the  characteristic  equation 
is  approximately 

g 

1  +  - 

s(s  -  X, 

and  the  root  locus  for  varying  My  (indicated  by  the  square  symbols  and 
solid  line)  is  as  shown  in  the  sketch  below: 


(s  -  Zw) 


)(s  -  Mg)(s  -  Zj 


=  0 


(42) 
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For  some  value  of  there  are  four  roots  (-  ,  a  +  its)  and  the  un¬ 

coupled  root  Zy,  which  are  shown  in  the  sketch  by  the  triangle  symbols. 

Also,  in  the  sketch,  is  a  pair  of  zeroes,  which  are  shown  by  the  circle 
symbol.  When  My  is  not  zero  the  characteristic  equation  is  approximately 


»0«.  (s2  -  Xu  s  ♦  ) 

(s  -  Zy)(s  +  X1)  ^s2  -  2o  s  +  (a2  +  ud2)J 


and  the  root  locus  for  varying  My  is  as  shown  by  the  dotted  line 
(My  >  0)  and  by  the  dashed  line  (My  <  0). 

When  the  wing  incidence  varies,  a  family  of  loci  results  as  shown  below: 
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The  locus  of  roots  of  the  unstable  mode  which  result  for  the  case  when  My 

equals  zero  during  a  transition  is  shown  by  a  double  dashed  line.  At  high 
incidence  angles,  the  actual  locus  indicates  a  somewhat  larger  instability 
than  the  My  equals  zero  case  due  to  the  small  positive  value  of  My. 

When  the  wing  incidence  is  less  than  roughly  hcP  ,  the  unstable  mode  becomes 
markedly  more  stable  because  of  the  rapid  increase  of  Mw  with  incidence 
in  a  statically  stable  sense  (negative),  as  shown  in  Figure  2k. 


From  this  brief  discussion,  it  can  be  seen  that  the  transition  from 
"helicopter-like"  motions  to  "airplane-like"  motions  is  strongly  influenced 
by  the  sign  and  magnitude  of  My.  For  example,  at  iw  equal  to  35  degrees, 

if  My  is  doubled,  the  dynamics  will  be  very  similar  to  those  of  a  normal 
airplane ;  and  if  at  iw  equal  to  20  degrees,  My  is  set  equal  to  zero, 

the  characteristic  modes  will  be  very  similar  to  a  conventional  helicopter 
(Figure  30). 


The  horizontal  center-of- gravity  location  will  influence  the  tern  Mw  as 

the  forward  speed  increases  and  the  variation  of  wing  lift  with  vertical 
velocity  becomes  significant.  In  the  design  stage,  then,  some  control  over 
this  derivative  is  available  at  higher  trim  speeds.  Various  programmed 
tail  incidence  settings  have  been  suggested  for  these  aircraft,  and  it 
should  be  realized  that  if  the  tail  is  programmed  in  such  a  way  that  it 
inadvertently  stalls  at  some  forward  speed  condition,  then  a  marked  un¬ 
favorable  change  in  My,  can  result  because  of  the  loss  of  the  horizontal 

tail  contribution  (Figure  2b).  It  may  be  difficult  to  estimate  when 
stalling  may  occur  due  to  lack  of  detailed  knowledge  of  downwash  angles  at 
the  tail.  The  possible  changes  in  this  derivative  due  to  stalling  of  the 
horizontal  tail  are  of  considerably  more  significance  than  any  effects  that 
tail  incidence  may  have  on  Mu  because  My  has  a  stronger  effect  on  the 

initial  stages  of  the  response. 


Figure  26  indicates  a  trend  of  the  characteristic  modes  that  would  be 
physically  expected  in  changing  from  hovering  to  forward  flight.  That  is, 
the  unstable  oscillation  in  hovering  becomes  the  conventional  phugoid  mode 
in  forward  flight,  and  the  two  convergences  present  in  hovering  form  the 
short-period  mode  in  forward  flight.  A  somewhat  more  peculiar  transition 
of  the  modes  occurs  if  the  lift  curve  slope  Zy  or  the  pitch  damping  Mg 

of  the  vehicle  remains  very  small,  as  may  be  seen  from  Figures  28  and  29. 
In  these  cases,  the  unstable  oscillation  becomes  the  short-period  mode  in 
forward  flight  and  the  two  real  roots  become  the  phugoid.  This  is  proba¬ 
bly  an  unusual  tendency  since  for  any  reasonable  vehicle  configuration  the 
lift  curve  slope  would  increase  with  forward  speed.  If  the  lift  curve 
slope  remains  zero,  the  airplane  would  certainly  have  peculiar  handling 
qualities  from  other  considerations.  Also,  if  the  angular  damping  is 
small,. the  dynamics  would  be  poor. 


To  summarize  the  results  of  this  section,  the  following  additional 
are  made. 
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With  respect  to  the  unstable  mode  that  becomes  the  conventional  phugoid 
mode  as  the  trim  speed  increases,  the  stability  derivative  is  very 

important  in  determining  the  motion  characteristics,  until  flight  con¬ 
ditions  are  encountered  where  My  becomes  large.  The  large  values  of 

My  (in  this  case  iw  s  2CP  )  tend  to  wash  out  the  effects  of  Mu.  The 
damping  of  this  mode  is  primarily  determined  by  Xy  and  M^.  The  angu¬ 
lar  damping  Mq  becomes  less  important  as  the  period  becomes  longer. 

At  low  speeds,  the  other  modes  consist  of  two  convergences;  the  large  root 
is  determined  by  Xy  and  Mq  and  the  small  root  by  Zy.  It  is  a  peculi¬ 
arity  of  these  aircraft  at  low  speeds  that  IxJ  >  I M^ | ,  which  is  opposite 

from  the  usual  case  of  a  helicopter  (Reference  12).  These  two  real  roots 
combine  to  become  an  oscillatory  motion  as  My  increases  with  trim  speed. 

This  oscillatory  mode  is  primarily  determined  by  My  ,  M^  and  Zy  ,  as 

would  be  expected  on  the  basis  of  the  conventional  short-period  approxi¬ 
mation  (Reference  11). 

COMPARISON  OF  THEORY  AND  EXPERIMENT 


Comparisons  of  the  theory  presented  here  and  experimental  data  measured  on 
a  0.10  scale  model  using  the  Princeton  Dynamic  Model  Track  are  shown  in 
Figures  32  through  36.  A  detailed  discussion  of  the  experimental  results 
appears  in  Reference  20.' ^'in  particular,  various  nonlinearities  appeared 
in  the  experimental  results  that  are  not  considered  here. 

The  comparison  of  trim  velocity  and  wing  incidence  is  shown  in  Figure  32. 
The  model  exhibited  higher  trim  velocities  than  calculated  by  the  theory. 
Inspection  of  Equation  (20b),  which  is  an  approximation  to  the  trim  ve- 
locity/wing  incidence  relationship,  indicates  either  that  the  flap  ef¬ 
fectiveness  used  in  the  calculations  was  higher  than  actually  present  on 
the  model  or  that  the  wing  lift  curve  slope  in  the  slipstream  is  smaller 
than  the  theoretical  value. 

However,  as  noted  in  the  previous  section,  the  horizontal  velocity  deriva¬ 
tives  are  not  particularly  sensitive  to  the  calculation  of  trim  velocity. 
If  the  horizontal  velocity  derivatives  are  reasonably  linear,  then  it  is 
implied  that  they  are  not  strongly  dependent  upon  trim  velocity.  There¬ 
fore,  no  attempt  was  made  to  improve  this  trim  calculation  before  pro¬ 
ceeding  with  the  derivative  comparisons.  It  is  felt  that  considerably 
more  detailed  aerodynamic  information  would  be  desirable  before  proceeding 
into  further  specific  investigations  of  this  nature;  that  is,  to  decide 
exactly  which  terms  are  in  error. 

For  the  derivative  comparisons,  only  the  force  and  moment  variations  with 
horizontal  velocity  are  considered  here  since  these  were  all  reasonably 
linear.  The  comparisons  are  shown  in  Figures  33,  31+  and  35.  While  the 
general  trend  of  the  derivatives  is  predicted,  the  agreement  is  not  par¬ 
ticularly  good.  The  tern  is  predicted  reasonably  well  in  hovering, 
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1 ut  experimentally  decreases  considerably  faster  than  the  theoretical 
value.  T!ie  term  Xy  is  underestimated  in  hover  and  agrees  reasonably  well 

at  lower  incidence  angles.  The  tern  Zy  is  predicted  reasonably  well  at 
incidence  angles  above  6C P  .  The  general  indications  are  that  at  high  incl- 


3L. 

dence  angles,  — -  is  predicted  reasonably  well,  but  that  the  prediction 

becomes  poorer  as  the  trim  speed  is  increased.  This  is  probably  due  to  the 
term  arising  from  the  lift  curve  slope  variation  which  becomes  quite  large 
at  low  incidence  as  well  an  increased  areas  of  separation  on  the  wing, 
arising  from  the  increase  in  average  effective  wing  angle  of  attack 


(Figure  5).  Thus,  it  appears  that  — —  is  overestimated.  It  is  also 

^u 

possible  that  the  contributions  of  the  horizontal  tail  to  the  derivative 
My  are  considerably  different  than  indicated  by  Figure  23.  That  is,  the 

slipstream  may  produce  dynamic  pressures  at  the  tail  that  are  considerably 
greater  than  free  stream,  causing  larger  effects  than  shown  in  Figure  23. 
The  experimental  values  of  My  presented  differ  somewhat  from  those 

obtained  in  a  preliminary  analysis  of  the  data  made  by  LTV.  Detailed 
analysis  of  the  static  and  dynamic  data  (Reference  20)  by  Princeton  indi¬ 
cates  that  My  is  negative  in  the  vicinity  of  6cP  wing  incidence  as 
shown. 


The  experimental  and  theoretical  comparison  of  the  pitch  damping  is  shown 
in  Figure  36.  The  theoretical  and  experimental  comparison  was  made  with¬ 
out  tail  rotor.  The  agreement  with  the  general  trend  of  this  derivative 
is  good.  While  at  very  low  speeds,  the  percentage  error  is  large,  the 
actual  derivative  is  so  small  that  it  is  probably  unreasonable  to  expect 
better  agreement,  and  the  correct  order  of  magnitude  is  obtained. 


No  comparison  is  made  of  the  derivatives  due  to  vertical  velocity  because 
of  the  nonlinearities  appearing  in  this  data.  The  majority  of  this  theo¬ 
retical  work  was  completed  prior  to  a  detailed  investigation  of  the  experi¬ 
mental  data  (Reference  20).  Further  remarks  on  the  comparison  of  this 
theory  with  experiment  will  be  found  in  Reference  20. 
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CONCLUSIONS 
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For  the  example  configuration  studied,  a  tilt-wing  VTOL  aircraft  of  the 

XC-142A  class  with  double  slotted  flaps,  and  the  assumptions  made,  the 

following  general  conclusions  are  stated: 

1.  In  order  to  predict  accurately  the  trim  velocity /wing  incidence  re¬ 
lationship  and  the  pitching  moment  required  to  trim  during  transitions 
of  tilt-wing  VTOL  type  vehicles,  it  is  necessary  to  have  precise 
knowledge  of: 

a.  The  aerodynamic  characteristics  of  the  wing/flap  combination 
with  the  authentic  slipstream/relative  wind  flow  field 
representation. 

b.  The  aerodynamic  characteristics  of  the  horizontal  tail,  with 
attention  to  the  downwash  behavior  and  dynamic  pressure 
variations . 

2.  Correlation  with  experimental  data  may  ;e  poor  because  of  the  assumed 
values  of  the  wing  lift  curve  slope  and  the  flap  effectiveness.  It 
appears  on  the  basis  of  the  predicted  "wing  incidence/forward  speed" 
relationship  and  the  importance  of  the  wing  contributions  to  the  sta¬ 
bility  derivatives  Xu  and  My  that  the  assumed  wing  lift  curve 
slope  is  too  high  in  the  wing  incidence  range  between  7CP  and  4cP  . 

Data  available  after  this  report  was  written  indicate  that  the  assumed 
flap  effectiveness  did  not  correspond  well  with  the  experimental  value. 

3.  On  the  basis  of  the  analysis,  the  propeller  inplane  force  and  wing 
drag  appear  to  have  relatively  small  contributions  to  the  stability 
derivatives . 

4.  The  slipstream/wing  lift  variations  dominate  the  changes  in  the  vehicle 
stability  derivatives  at  slow  speeds.  The  propeller  thrust  force  vari¬ 
ations  do  give  appreciable  contributions,  but  they  appear  to-be  rela¬ 
tively  less  significant  when  compared  to  the  gross  changes  in  the  de¬ 
rivatives  caused  by  the  slipstream/wing  lift  variations.  TKi!£;'par- 
ticular  relative  proportion  of  changes  is  due  to  some  extent  to  the 
layout  design  or  geometry  of  the  vehicle  studied,  but  the  example  con¬ 
figuration  geometry  is  generally  typical  of  tilt-wing-type  aircraft. 

5.  On  the  basis  of  a  limited  comparison  with  experimental  results  from  a 
scale  model  of  the  XC-142A,  the  important  stability  derivatives  were 
predicted  quite  well  in  hovering,  and  the  general  trend  with  decreasing 
incidence  was  predicted;  however,  the  magnitudes  predicted  are  con¬ 
siderably  in  error  for  wing  incidences  representing  the  region  of 
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flight  midway  through  transition  (around  45  degrees).  The  assumed 
configuration  is  somewhat  dissimilar  from  the  scale  model. 

6.  Favorable  effects  on  t.ne  attitude  stability  (My)  and  pitch  damping 
(M^)  arise  from  the  increasing  dynamic  pressure  at  the  tail  as  wing 
incidence  is  reduced. 

7.  The  change  in  stability  characteristics  from  "helicopter-like"  dynamics 
at  high  incidence  to  "airplane-like"  dynamics  at  low  incidence  is 
largely  determined  by  the  magnitude  and  sign  of  the  attitude  stability 
(My)  of  the  vehicle. 
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RECOMMENDATIONS 


1.  Further  comparison  of  this  theory  with  experiment  is  desirable.  It 
would  be  valuable  to  pursue  comparison  of  the  expressions  developed 
with  various  experimental  data  on  isolated  components  available  be¬ 
fore  proceeding  with  further  comparisons  of  complete  configurations. 

2.  Due  to  the  increasing  importance  of  the  horizontal  tail  in  providing 
favorable  effects  as  the  incidence  decreases,  experimental  information 
on  downwash  characteristics  in  the  vicinity  of  the  horizontal  tail  of 
tilt-wing  configurations  is  desirable. 

3.  Consideration  should  be  given  to  the  importance  of  nonlinearities  in 
regions  where  the  linearized  theory  predie4-..  derivatives  near  zero, 
as  for  example  in  the  region  where  the  alt  .de  stability  (M,^) 
changes  sign. 

U.  Since  there  is  at  present  considerable  experimental  information  on 

the  aerodynamic  characteristics  of  wing -propeller  combinations,  these 
data  should  be  organized  in  such  a  fashion  as  to  verify  the  semi- 
empirical  approach  used  to  predict  wing  lift  curve  slope,  or  to  suggest 
modification  to  the  form  used  here. 
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STEADY  LEVEL  FLIGHT ,  a  *0* 


FIGURE  I.  TRIM  VELOCITY  VS.  WING  INCIDENCE  ANGLE  FOR 
VARIOUS  VALUES  OF  FLAP  EFFECTIVENESS. 


INDUCED  VELOCITY  DIVIDED  BY  TIP  SPEEO  (X,) 


FIGURE  2.  PROPELLER  THRUST  AND  WING  LIFT  VS.  WING 
INCIDENCE. 
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FIGURE  3.  PROPELLER  BLADE  ANGLE  AND  NON  DIMENSIONAL 
INDUCED  VELOCITY  VS.  WING  INCIDENCE. 
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figure  4.  SLIPSTREAM  VELOCITY  VS.  WING  INCIDENCE  FOR 
VARIOUS  VALUES  OF  FLAP  EFFECTIVENESS. 


STEADY  LEVEL  FLIGHT,  a  -  O* 


FIGURE  5.  EFFECTIVE  WING  ANGLE  OF  ATTACK  VS.  WING 
INCIDENCE. 
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FIGURE  6.  PITCHING  moment  required  to  trim 
VS.  WING  incidence. 
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FIGURE  7.  INDUCED  VELOCITY  DERIVATIVES  VS.  WING 
INCIDENCE. 
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STEADY  FLIGHT,  Q  "0* 


FIGURE  9.  DEPENDANCE  OF  VARIOUS  SUPSTREAM  PARAMETERS  ON 
WING  NCIDENCE  /  TRIM  VELOCITY  RELATIONSHIP. 
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*  Gu*f  -0  RATE  OF  CHANGE  OF  SLIPSTREAM  VELOCITY  WITH 
'OPi/ONTAl  VELOCITY  vs  WING  'NClDENCE 


»iGjnr  rate  of  change  of  effective  wing  angle  of 

ATTACK  with  HORIZONTAL  VELOCITY  *S  WING  1NC1OENCE 


r<au*C  2  RATE  OF  CHANGE  OF  TmRUST  WIth  HORIZONTAL 
VELOCITY  VS  WING  INCIDENCE" 
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RE  13. 


RATE  OF  CHANGE  OF  WING  LIFT  WITH 
HORIZONTAL  VELOCITY  VS.  WING  INCIDENCE. 


RATE  OF  CHANGE  OF  HORIZONTAL  FORCE  WITH 
HORIZONTAL  VELOCITY  DIVIDED  BY  MASS;  PER  SECOND 


WING  INCIDENCE;  DEGREES  (i*) 


FIGURE  14.  STABILITY  DERIVATIVES  -  RATE  OF  CHANGE  OF 
HORIZONTAL  FORCE  WITH  HORIZONTAL  VELOCITY 
VS.  WING  INCIDENCE. 


44 


FIGURE  15.  STABILITY  DERIVATIVES  -  RATE  OF  CHANGE  OF 
VERTICAL  FORCE  WITH  HORIZONTAL  VELOCITY 
VS.  WING  INCIDENCE. 


FtQU«  *•  RATE  OF  CHANGE  OF  PROPELLER  THRUST  WITH 
VERTICAL  VELOCITY  VS  WING  INCIDENCE 
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figure  20.  STABILITY  DERIVATIVES  -  RATE  OF  CHANGE  OF 

horizontal  force  with  vertical  VELOCITY 

VS  WING  INCIDENCE. 
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FIGURE  21.  STABILITY  DERIVATIVES  -  RATE  OF  CHANGE  OF 
VERTICAL  FORCE  WITH  VERTICAL  VELOCITY 
VS.  WING  INCIDENCE. 
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STEADY  LEVEL  FLIGHT,  C!  *0* 
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FIGURE  22.  STABILITY  DERIVATIVES  -  RATE  OF  CHANGE  OF 
PITCHING  MOMENT  WITH  HORIZONTAL  VELOCITY 
VS  WING  INCIDENCE  .  ( NO  HORIZONTAL  TAIL  ) 
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RATE  OF  CHANGE  OF  PITCHING  MOMENT  WITH  HORIZONTAL 
VELOCITY  DIVIDED  BY  INERTIA;  PER  SECOND  SQUARED 
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(EQ.  39,  p.23) 


FIGURE  23.  HORIZONTAL  TAIL  CONTRIBUTIONS  TO  PITCHING 
MOMENT  CHANGE  WITH  HORIZONTAL  VELOCITY. 
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FIGURE  24a.  STABILITY  DERIVATIVES  -  RATE  OF  CHANGE  OF 
PITCHING  MOMENT  WITH  VERTICAL  VELOCITY 
VS.  WING  INCIDENCE. 
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FIGURE  27.  LOCUS  OF  ROOTS  AS  A  FUNCTION  OF  WING  INCIDENCE 
SENSITIVITY  TO  Xy  • 


FIGURE  29.  LOCUS  OF  ROOTS  AS  A  FUNCTION  OF  WING  INCIDENCE 
SENSITIVITY  TO  M#  • 


FIGURE  30.  LOCUS  OF  ROOTS  AS  A  FUNCTION  OF  WING  INCIDENCE 
SENSITIVITY  TO  M*. 


FIGURE  31.  LOCUS  OF  ROOTS  AS 
SENSITIVITY  TO  M„  . 
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FIGURE  32. 


""N 


STEADY  LEVEL  FLIGHT, a  •  <>• 


WING  INCIDENCE!  DEGREES 


(i») 


COMPARISON  OF  THEORY  AND  EXPERIMENT 
TRIM  VELOCITY  VS.  WING  INCIDENCE. 


STEADY  LEVEL  FLIGHT,  a  -0* 


WING  INCIDENCE;  DEGREES  (iw) 


FIGURE  33.  COMPARISON  OF  THEORY  AND  EXPERIMENT  - 

Xg  VS.  WING  INCIDENCE. 
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FIGURE  34  COMPARISON  OF  THEORY  AND  EXPERIMENT  - 

Zu  V&  WING  INCIDENCE. 
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FIGURE  35.  COMPARISON  OF  THEORY  AND  EXPERIMENT - 

Mu  VS.  WING  INCIDENCE. 
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FIGURE  36.  COMPARISON  OF  THEORY  AND  EXPERIMENT  - 
Mtf  VS.  WING  INCIDENCE 
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FIGURE  37.  ASSUMED  FLAP  CHARACTERISTICS. 
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APPENDIX  A 


PHYSICAL  PARAMETERS  OF  ASSUMED  AIRCRAFT  AND 
"DESCRIPTION  OF  APPARATUS  AND  EXPERIMENTS 


ASSUMED  AIRCRAFT 

W  =  37,350  pounds 
Propellers 
N  =  4 

R  =  7-75  feet 

fj  =  129  radians  per  second 

Wing 

S  a  534  square  feet 
bw  =  67.5  feet 
c  =  8.07  feet 
AR  =  8.53 
K  =  0.5 

Horizontal  Tail 

3^  =  140  square  feet 

AR  =  5.08 


Tail  Rotor 

R1  =  4  feet 

0  =  363*4  radians  per  second 

T 

a  =  0.124 
T 


Iy  =  100,000  slug -feet  squared 

ct  =  0.153 
b  =  4 

a  =  5>73  per  radian 
Cj^  =  5*0  per  radian 

C  =  0.0130 

uo 

qtq  =  0.027  radian 

s6 

— £  =  0.333  (double-slotted  flap) 

s 

l  =  24.3  feet 
HT 

t  =  6.16  feet 

H 

^TR  =  31.4  feet 

b  =  3 

a  =  5*73  per  radian 

T 
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The  flap  effectiveness  is  shown  in  Figure  37-  The  angle  of  zero  lift  of 
the  wing  (crQ),  the  zero  lift  drag  coefficient  (C_  ),  and  the  pitching 

moment  about  the  aerodynamic  center  (C^)  are  assumed  to  vary  as  shown 

with  flap  deflection.  The  flap  deflection  is  programmed  as  a  function  of 
wing  incidence  (Figure  A-2). 

Various  geometric  distances  of  the  aircraft,  ,  -ty  ,  -t-oP  >  >  which 

are  measured  from  the  center  of  gravity  of  the  airplane,  as  shown  in 
Figure  A-lb,  are  expressed  analytically  (in  feet)  ar  follows; 


0.158  +  (0.19)  sin  iw  +  (2.61)  cos  iw 
1.105  +  (2.61)  sin  iv  +  (0.19)  cos  iy 
1.06  +  l 

OW 


These  quantities  are  shown  graphically  in  Figure  A-3. 

DESCRIPTION  OF  APPARATUS  AND  EXPERIMENTS 

The  experimental  results  presented  were  obtained  from  a  0.10  scale  dynamic 
model  of  the  XC-112A  tested  on  the  Princeton  Dynamic  Model  Track. 

Test  Facility 


lhe  Princeton  University  Dynamic  Model  Track  is  a  unique  facility,  designed 
expressly  for  the  study  of  the  dynamic  motions  of  helicopter  and  VTOL 
models  at  equivalent  flight  speeds  of  up  to  60  knots  (for  a  one-tenth  scale 
model).  Basic  components  of  the  facility  include  a  servo-driven  carriage 
riding  on  a  track  750  feet  long,  located  in  a  building  of  cross  section 
30  x  30  feet;  the  carriage  has  an  acceleration  potential  of  0.6g  and  a 
maximum  speed  of  10  feet  per  second.  A  detailed  description  of  the  facili¬ 
ty  and  the  testing  techniques  employed  may  be  found  in  Reference  22. 

A  model  may  be  attached  to  the  carriage  by  one  of  several  booms.  One  mount 
permits  relative  displacements  of  the  model  with  respect  to  the  carriage  in 
horizontal  and  vertical  directions,  as  well  as  allowing  it  to  rotate  in  the 
plane  determined  by  these  two  directions.  Horizontal  relative  motion  of 
the  model  with  respect  to  the  carriage  is  sensed  and  used  to  command  the 
carriage  to  follow  the  model  in  a  closed  loop  fashion.  Similarly,  vertical 
displacement  of  the  model  with  respect  to  the  carriage  commands  the  boom 
to  move  vertically.  This  servo  operation  of  the  carriage  allows  the  model 
to  fly  "free"  with  no  restraints  on  the  dynamic  motions  being  investigated. 


70 


t 


...  1 


This  method  of  testing  may  be  considered  similar  to  dynamic  flight  testing, 
with  considerably  more  control  over  the  experiment  possible. 

In  addition  to  the  dynamic  testing  as  described  above,  atic  testing  is 
conducted  by  programming  carriage  movement  in  accordance  with  pre-selected 
velocity  profiles  at  constant  angle  of  attack.  Programmed  angle-of-attack 
changes  at  constant  velocity  are  also  possible.  The  model  is  rigidly 
mounted  on  v.he  carriage  and  forces  and  moments  acting  on  the  model  are 
measured  with  strain  gauges.  Although  these  static  experiments  are  similar 
to  wind  tunnej.  tests,  this  facility  offers  the  advantages  of  a  30  x  30 
foot  test  section  with  a  uniform  air  velocity,  free  from  turbulence.  Pre¬ 
cise  speed  control  over  a  range  of  speeds  from  backward  flight  through 
hover  to  forward  flight  is  available. 

Model 

A  three-view  drawing  of  the  0.10  scale  dynamic  model  of  the  XC-142A  is 
shown  in  Figure  A-4. 

The  fuselage  is  constructed  of  an  inner  and  outer  Fiberglas  skin,  vacuum 
molded  and  bonded  to  a  styrofoam  core.  An  aluminum  box  spar  is  the  main 
structural  member  of  the  wing.  Mahogany  ribs  and  a  vacuum  molded  Fiberglas 
wing  surface  form  the  external  airfoil  shape.  Tne  double-slotted  flaps 
were  constructed  of  low  density  styrofoam  with  a  Fiberglas  covering. 

The  model  drive  motor  is  a  200-volt,  400-cycle,  three-phase  electric  motor, 
rated  at  five  horsepower,  mounted  on  a  bulkhead  in  the  fuselage.  Power  for 
the  four  propellers  was  transmitted  to  a  central  transmission  and  from 
thence  to  right-angle  gear  boxes  located  in  the  wing  by  flexible  shafting. 

A  separate  power  take-off  is  used  to  drive  the  tail  rotor.  Propeller  gear 
boxes  and  housings  were  mounted  directly  on  the  wing  spar.  The  propeller 
blades  were  constructed  of  Fiberglas  by  the  Hamilton  Standard  Division  of 
;he  United  Aircraft  Corporation. 

Model  control  positions  are  set  from  a  control  console  on  the  carriage. 

The  model  incorporated  electrically  controllable  blade  angles  on  each  of 
the  four  propellers.  The  blade  angle  of  the  tail  rotor  is  also  variable 
to  provide  pitching  moment  trim.  Wing  incidence,  flaps,  ailerons,  and  the 
horizontal  tail  are  also  power  operated  so  that  transition  runs  may  be 
made  with  selected  programming  of  all  required  controls.  All  of  these 
systems  are  closed  loop  position  controls.  The  response  characteristics 
of  each  of  these  loops  is  sufficiently  fast  to  permit  study  of  automatic 
stabilization  systems. 

Static  Data 


Measurements  of  vertical  force  (perpendicular  to  the  free  stream),  hori¬ 
zontal  force  (parallel  to  the  free  stream),  and  pitching  moment  acting  on 
the  model  near  trim  conditions  (X  =  0)  were  made  to  determine  the  static 
stability  derivatives  of  this  vehicle  at  various  wing  incidences.  The  flap 
deflection  is  given  as  a  function  of  wing  incidence  in  Figure  A-2. 
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Fuselage  angle  of  attack,  propeller  blade  angle,  and  propeller  RPM  were 
held  constant,  and  the  carriage  was  programmed  for  a  very  small  acceler¬ 
ation  (*  O.Olg's)  such  that  the  velocity  of  the  model  was  varied  about  a 
trim  condition  (X  >0).  This  velocity  program  was  conducted  at  three 
fuselage  angles  of  attack  (a«0,  a  as  ♦  15°,  o  a  -  15° ) ,  one  propeller 
rotational  speed  (4000  RPM),  and  one  blade  angle  (17.5°)  setting.  At  a 
wing  incidence  of  40*,  some  data  were  taken  at  two  other  blade  angle 
settings. 

Previous  experience  with  this  technique  of  quasi-static  testing  has  shown 
that  the  derivatives  are  identical  to  those  obtained  from  point  by  point 
measurements  at  velocity  increments  about  a  trim  condition.  A  considerable 
reduction  in  testing  time  is  realized  provided  that  the  carriage  acceler¬ 
ations  involvsd  are  kept  very  small. 

The  values  of  the  stability  derivatives  determined  from  these  experiments 
are  given  in  Figures  32  through  36.  Detailed  consideration  of  these  data 
may  be  found  in  Reference  20. 


b.  GEOMETRIC  DISTANCES  FROM  CENTER  OF  GRAVITY 


—  Proptlltf  Ptf  n* 
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APPENDIX  B 

ANALYTICAL  EXPRESSIONS  FOR  STABILITY  DERIVATIVES 


EqJATIONS  OF  MOTION 


The  equations  of  motion  with  respect  to  stability  axes  for  motions  in  the 
plane  of  symmetry  of  an  aircraft  (Reference  11)  may  be  written  as  follows: 

m(u  +  w  0)  =  X(u,w)  -  W  sin  0  (B-l) 

ra(w  -  u  0)  =  Z(u,w)  -  W  cos  @  (B-2) 

I  0  =  M(u,w,0)  (B-3) 


Expanding  the  aerodynamic  forces  X,  Z,  and  M  in  a  Taylor  series  about 
level,  unaccelerated  flight,  retaining  only  first-order  terms,  and 
linearizing  the  inertia  and  gravity  terms  about  0O  =0,  uQ  =  VQ,  wQ  =  0, 
result  in  the  following: 


(B-4) 

(B-5) 

(B-6) 


The  stability  derivatives  have  been  calculated  using  the  following  ex¬ 
pressions,  based  on  the  assumptions  discussed  in  the  text. 


STABILITY  DERIVATIVES 


Resolving  the  propeller  and  wing  forces  in  the  horizontal  direction  as 
shown  in  Figure  A-l,  the  X-force  may  be  written  as 


X(u,w)  =  T  cos  iw  -  H  sin  iw  -  Ls  sin(iw  -  ae )  -  Ds  cos(iw  -  <*e)  +  (x)f 

(B-7) 
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Taking  derivatives  of  the  above  expression,  the  horizontal  force  deriva¬ 
tives  are 
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Resolving  farces  in  the  vertical  direction,  the  Z-force  is 


Z(u,w)  =  -  T  sin  iy  -  H  cos  i¥  -  Lg  coe(iw  -  or^.)  +  sin(iw  -  ag)  (B-JO) 


and  the  Z-force  derivatives  are 
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The  pitching  moment  acting  on  the 
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may  be  expressed  as 
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The  pitching  moment  derivatives  are  obtained  by  taking  derivatives  of 
this  expression: 
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Eq.  (B  -15)  cont. 
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PROPELLER  EQUATIONS 
Induced  Velocity 

The  average  value  of  the  propeller  induced  velocity  may  be  expressed  in 
terms  of  the  propeller  blade  angle  and  the  velocities  at  the  propeller 
plane  by  eliminating  the  thrust  coefficient  from  the  momentum  equation 
using  the  blade  element  equation  for  thrust  coefficient. 
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This  equation  may  be  rearranged  to  give  a  fourth  order  equation  for  X 
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The  derivatives  of  the  induced  velocity  are  obtained  from  Equation  (B-18) 
as 


(B-19) 


(B-20) 


Thrust  Equations 

The  blade  element  expression  for  the  thrust  (Reference  10,  page  190)  is 

T  .  p  n  R2  (OR)2  22  r6-^  (1  +  as! )  .  h_l!l  (B-21) 

2  L  3  2  2  J 


This  expression  is  written  in  terms  of  nondimensional  velocities  perpen¬ 
dicular  and  parallel  to  the  propeller  plane  (v,n).  To  convert  to  the 
velocities  u  and  w  used  in  the  stability  analysis,  the  following  re¬ 
lationships  are  used: 
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and 


w  =  V  sin  a 


u  =  V  cos  a 


(B-23) 


Therefore,  the  relationships  between  derivatives  are 
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Pitching  Moment  Equations 

The  propeller  pitching  moment,  arising  from  the  variation  in  induced  ve¬ 
locity  across  the  disc,  is 
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where,  from  Reference  8, 
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and  the  pitching  moment  variations  are 
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where,  from  Equation  (B-32), 
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and  from  Equations  (B-2U)  and  (B-24a),  and  using  the  definition  of  X^, 
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The  pitching  moment  arising  from  an  angular  velocity  about  the  hub  is 


<m6>id  =  -  pg"nB5 


(B-39) 


when  the  induced  velocity  has  been  assumed  constant,  as  discussed  in 
the  text . 

WUffG  EQUATIONS 

"  * 

The  wing  is  assumed  to  be  completely  immersed  in  the  slipstream,  and  so 
only  those  contributions  are  given. 
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Slipstream  Velocity 


With  the  assumption  that  the  wing  forces  are  dependent  upon  the  fully 
developed  induced  velocity,  the  slipstream  velocity  is  expressed  as 
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The  slipstream  velocity  derivatives  are 
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Effective  Wing  Angle  of  Attacx 

The  effective  wing  angle  of  attack  in  the  slipstream  may  be  expressed  as 
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The  effective  angle  of  attack  derivatives  are 
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where  the  slipstream  velocity  derivatives  are  given  by  Equations  (B-4l) 
and  (B-42). 

Wing  Lift 

The  wing  lift  in  the  slipstream  is  expressed  as 

Ls  -  I  P  S  CL,S  VR  (°o  +  «e>  (E-U5) 


where 


c  =  c  r  C 

L.S  L.S  v  L 

'  or  ’  a  a 


i 

w 


or 


) 
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where  the  effective  angle  of  attack  and  slipstream  velocity  derivatives 
are  given  above.  The  lift  curve  slope  derivatives  are  listed  below. 

Wing  Drag 


The  wing  drag  in  the  slipstream  is  expressed  as: 
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The  derivatives  are 
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Wing  Lift  Curve  Slope  in  Slipstream 

As  discussed  in  the  text,  it  is  assumed  that  the  wing  lift  curve  slope  in 
the  slipstream  may  be  expressed  aE 
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TAIL  ROTOR  EQUATIONS 


The  calculation  of  the  tail  rotor  forces  is  based  on  the  same  assumptions 
as  used  for  the  propellers. 
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The  nondimensional  velocities  perpendicular  and  parallel  to  the  plane  of 
the  tail  rotor  (taken  parallel  to  the  fuselage  reference  here)  are 
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The  downwash  angle,  e^,  ,  arises  from  the  lift  of  the  horizontal  tail  and 
is  expressed  as 
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The  relationship  between  the  derivatives  with  respect  to  the  yelocities 
used  in  the  stability  analysis  and  the  derivatives  taken  with  respect  to 
velocities  perpendicular  and  parallel  to  the  plane  of  the  tail  rotor  is 
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where,  from  Equation  (A-5*0> 
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5XlT  and  can  be  calculated  by  the  same  method  as  that  of 
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The  transformation  equations  are 
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The  inplane  force  of  the  tail  rotor  and  the  pitching  moment  of  the  tail 
rotor  about  its  hub  are  neglected. 

HORIZONTAL  TAIL  EQUATIONS 


The  horizontal  tail  lift  is  expressed  as 
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The  derivatives  of  this  expression  are 
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The  drag  of  the  horizontal  tail  is  neglected. 


(B-72) 


(B-73) 


(B-74) 


90 


V' 


Unclassified _ 

Security  Cla»si(ication 


DOCUMENT  CONTROL  DATA  •  RAD 


(Security  claaat/lcation  ot  tltla.  body  of  abstract  and  mdauin§  annotation  muat  b  a  antarad  wh  an  tha  ova  rail  raport  la  c  laaailtad) 

1  ORIGINATING  ACTIVITY  (Corporata  author)  2a  RCRORT  IICuRlTv  c  LAIHRiCA  TiON 

'npT^n'r+.mMrrf-.  n-T  AprnQnnpp  nnrl  Mpnhftn 1 pn 1  pnnps  Unclassified 


Department  of  Aerospace  and  Mechanical  Sciences 
Princeton  University 


26  QROUF 


1  REPORT  TITLE 

An  Analytical  Study  of  Factors  Influencing  the  Longitudinal  Stability  of 
Tilt-Wing  VTOL  Aircraft 

4  DESCRIPT  IVE  NOrES  (Typa  ot  raport  and  Inctu&i  va  a  *taa) 

5  AUTHORfS)  (La*f  nama.  tirat  nama.  Initial) 

Eeppu,  G.  and  Curtiss,  H.  C.,  Jr. 


e  RE PO  RT  DATE 

July  1906 _ 

8a  CONTRACT  OR  GRANT  NO 

DA  UL-177-AMC-8(T) 

b.  PROJECT  NO 

1P125901A14233 


10  AVA  (LABILITY /  LIMITATION  NOTICES 


7a  TOTAL  NO.  OR  PAGCf 


7b  NO  OF  REF* 


9a  ORIGINATOR’S  REPORT  NUMB E A(S) 


USAAVLABS  Technical  Report  66-53 


9b  OTHER  REPORT  NOfSj  (Any  othar  numbara  that  may  ba  aaal&ad 
ihla  raport ) 

Princeton  Uniy.  Aero.  &  Mech.  Sci. 

Report  No.  756 


Distribution  of  this  document  is  unlimited. 


It.  SUPPL  EMENTARY  NOTES 


I)  ABSTRACT 


12  SPONSORING  MILITARY  ACTIVITY 


U.  S.  Army  Aviation  Materiel  Laboratories! 
Fort  Eustis,  Virginia 


An  analytical  method  for  predicting  the  stability  characteristics  of  tilt-wing 
VTOL  aircraft  in  the  transition  speed  range  is  presented.  Sample  calculations 
based  on  an  assumed  tilt-wing  VTOL  transport  configuration  of  the  XC-142A  class 
with  double -slotted  flaps  are  given. 

Particular  emphasis  is  placed  on  the  sensitivity  of  the  results  to  various 
assumptions  made  in  the  analysis.  The  contributions  of  the  various  aircraft 
components  and  the  aerodynamic  interactions  of  the  components  to  the  stability 
derivatives  are  discussed,  as  well  as  the  changes  in  the  characteristic  modes 
of  motion  of  the  vehicle  that  result  from  variations  in  the  stability  derivatives. 
The  trim  conditions  of  the  vehicle  are  shown  to  be  quite  sensitive  to  the 
prediction  of  the  flap  characteristics. 

A  limited  comparison  of  the  calculated  results  with  experimental  data  obtained 
from  a  dynamic  model  of  the  XC-142A,  which  is  somewhat  dissimilar  from  the 
assumed  configuration,  is  presented.  This  comparison  indicates  that  the  trends 
of  the  stability  derivatives  are  correctly  predicted.  The  agreement  between 
theory  and  experiment  is  good  in  hovering;  however,  as  the  wing  incidence  is 
reduced,  the  difference  between  theory  and  experir ent  becomes  quite  large. 


Unclassified 


Security  Classification 


key  aonos 


stability,,  longitudinal 
airplanes,  VTOL  -  stability 
airplanes,  tilt-wing 
airplanes  (XC-142A) 


INSTRUCTIONS 


I.  ORIGINATING  ACTIVITY:  Enter  the  name  and  address 
of  the  contractor,  subcontractor,  grantee,  Department  of  De¬ 
fense  activity  or  other  organization  (corporate  author)  issuing 
I  he  report. 

2a.  REPORT  SECUHTY  CLASSIFICATION:  Enter  the  over¬ 
all  security  classification  of  the  report.  Indicate  whether 
"Restricted  Data"  is  included.  Marking  is  to  be  in  accord¬ 
ance  with  appropriate  security  regulations. 

2h.  GROUP:  Automatic  downgrading  is  specified  in  DoD  Di¬ 
rective  5 200. 10  and  Armed  Forces  Industrial  Manual.  Enter 
the  group  number.  Also,  when  applicable,  show  that  optional 
markings  have  been  used  for  Group  3  and  Group  4  as  author¬ 
ized. 

3.  REPORT  TITLE:  Enter  the  complete  report  title  in  all 
capital  letters.  Titles  in  all  cases  should  be  unclassified. 

If  a  meaningful  title  cannot  be  selected  without  classifica¬ 
tion,  show  title  classification  in  all  capitals  in  parenthesis 
immediately  following  the  title. 

4.  DESCRIPTIVE  NOTES:  If  appropriate,  enter  the  type  of 
report,  e.g.,  interim,  progress,  summary,  annual,  or  final. 

Give  the  inclusive  dates  when  a  specific  reporting  period  is 
covered. 

5.  AUTHOR(S):  Enter  the  namefs)  of  authors)  as  shown  on 
or  in  the  report.  Enter  last  name,  first  name,  middle  initial. 

If  military,  show  rank  and  branch  of  service.  The  name  of 
the  principal  author  is  an  absolute  minimum  requirement. 

6.  REPORT  DATE:  Enter  the  date  of  the  report  as  day, 
month,  year:  or  month,  year.  If  more  than  one  date  appears 
on  the  report,  use  date  of  publication. 

7a.  TOTAL  NUMBER  OF  PAGES:  The  total  page  count 
should  follow  normal  pagination  procedures,  i. e. ,  enter  the 
number  of  pages  containing  information 

7b.  NUMBER  OF  REFERENCES  Enter  the  total  number  of 
references  cited  in  the  report. 

8a.  CONTRACT  OR  GRANT  NUMBER:  If  appropriate,  enter 
the  applicable  number  of  the  contract  or  grant  under  which 
the  report  was  written. 

8b,  8c,  ba  id  PROJECT  NUMBER:  Enter  the  appropriate 
military  department  identification,  such  as  project  number, 
subproject  number,  system  numbers,  task  number,  etc. 

9 a  ORIGINATOR'S  REPORT  NUMBER(S):  Enter  the  offi¬ 
cial  report  number  by  which  the  document  will  be  identified 
and  controlled  by  the  originating  activity.  Thir  -lumber  must 
be  unique  to  this  report. 

9b  OTHER  REPORT  NUMBER(S):  If  the  report  has  been 
assigned  any  other  report  numbers  (either  by  the  originator 
or  by  the  sponsor),  also  enter  this  numberfs). 


10.  AVAILABILITY/LIMITATION  NOTICES:  Enter  any  lim¬ 
itations  on  further  dissemination  of  the  report,  other  than  those 
imposed  by  security  classification,  using  standard  statements 
such  as: 

(1)  "Qualified  requesters  may  obtain  copies  of  this 
report  from  DDC.” 

(2)  "Foreign  announcement  and  dissemination  of  this 
report  by  DDC  is  not  authorized.” 

(3)  ”U.  S.  Government  agencies  may  obtain  copies  of 
this  report  directly  from  DDC.  Other  qualified  DDC 
users  shall  request  through 


(4)  "U.  S.  military  agencies  may  obtain  copies  of  this 

report  directly  from  DDC  Other  qualified  users 
shall  request  through 


(5)  "All  distribution  of  this  report  is  controlled.  Qual¬ 
ified  DDC  users  shall  request  through 


if  the  report  has  been  furnished  to  the  Office  of  Technical 
Services,  Department  of  Commerce,  for  sale  to  the  public,  indi¬ 
cate  this  fact  and  enter  the  price,  if  known. 

11.  SUPPLEMENTARY  NOTES:  Use  for  additional  explana¬ 
tory  notes. 

12.  SPONSORING  MILITARY  ACTIVITY:  Enter  the  name  of 
the  departmental  project  office  or  laboratory  sponsoring  (pay¬ 
ing  lor)  the  research  and  development.  Include  address, 

13  ABSTRACT  Enter  an  abstract  giving  a  brief  and  factual 
summary  of  the  document  indicative  of  the  report,  even  though 
it  mav  also  appear  elsewhere  in  the  body  of  the  technical  re¬ 
port  If  additional  space  is  required,  a  continuation  sheet 
shall  be  attached. 

Il  is  highly  desirable  that  the  abstract  of  classified  re¬ 
ports  be  unclassified.  Each  paragraph  of  the  abstract  shall 
end  with  an  indication  of  the  military  security  classification 
of  the  information  in  the  paragraph,  represented  as  (TS),  (S), 

(C).  or  (V). 

There  is  no  limitation  on  the  length  of  the  abstract.  How¬ 
ever,  the  suggested  length  is  from  150  to  225  words. 

14  KEY  WORDS:  Key  words  are  technically  meaningful  terms 
or  short  phrases  that  characterize  a  report  and  may  be  used  as 
index  entries  for  cataloging  the  report.  Key  words  must  be 
selected  so  that  no  security  classification  is  required.  Iden- 
fiers,  such  as  equipment  model  designation,  trade  name,  mili¬ 
tary  project  code  name,  geographic  location,  may  be  used  as 
key  words  but  will  be  followed  bv  an  indication  of  technical 
context.  The  assignment  of  links,  r  iles,  and  weights  is 
optional 


Unclassified 


